International Journal of Modern Physics A 
© World Scientific Publishing Company 



STRING QUANTUM GRAVITY, LORENTZ-INVARIANCE 
VIOLATION AND GAMMA-RAY ASTRONOMY* 



NICK E. MAVROMATOS 
CERN, Theory Division, CH-1211 Geneva 23, Switzerland; 
On leave from: Department of Physics, Theoretical Physics, King's College London, Strand 

London, WC2R 2LS, UK 



Received 17th August 2010 

In the first part of the review, I discuss ways of obtaining Lorentz-Invariance- Violating 
(LIV) space-time foam in the modern context of string theory, involving brane world 
scenarios. The foamy structures are provided by lower-dimensional background brane 
defects in a D3-brane Universe, whose density is a free parameter to be constrained phe- 
nomenologically. Such constraining can be provided by high energy gamma-ray photon 
tests, including ultra-high energy /infrared photon-photon scattering. In the second part, 
I analyze the currently available data from MAGIC and FERMI Telescopes on delayed 
cosmic photon arrivals in this context. It is understood of course that conventional As- 
trophysics source effects, which currently are far from being understood, might be the 
dominant reason for the observed delayed arrivals. I also discuss how the stringent con- 
straints from studies of synchrotron-radiation from distant Nebulae, absence of cosmic 
birefringence and non observation of ultra-high-cncrgy cosmic photons can be accommo- 
dated within the aforementioned stringy space-time foam model. I argue that, at least 
within the currently available sets of astrophysical data, the stringy foam model can 
avoid all theses constraints in a natural range of the string coupling and mass scale. The 
key features are: (i) transparency of the foam to electrons and charged probes in general, 
(ii) absence of birefringence effects and (iii) a breakdown of the local effective lagrangian 
formalism. However, in order to accommodate, in this theoretical framework, the data 
of the FERMI satellite on the delayed arrival of photons from the short intense Gamma 
Ray Burst GRB 090510, in a way consistent with the findings of the MAGIC telescope, 
a non uniform density of brane foam defects must be invoked. 
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1. Introduction and Summary 

An alternative title for the review article could be MA GIC Strings : Usually the ter- 
minology M-theory, with M standing for either Magical or Marvelous or Mysterious, 
is attributed to the underlying (yet not completely understood) unifying theory of 
all known string theories as a result of the many appealing duality and other 
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symmetries it possesses, which result in the unification of the five known string the- 
ories, viewed as a low-energy limit of M-theory. This is a super-unification picture, 
which may prompt the way for a detailed understanding of the yet elusive theory 
of the quantum structure of space-time, otherwise termed as "Quantum Gravity" 
(QG). Our current knowledge/understanding of M-theory is limited. Schematically, 
this knowledge is restricted to a few small regions of the interior of the diagram 
of fig. 1. Nevertheless, for an analysis of some of the predictions of string theory 
that could have some relevance to observable low-energy physics this may not be 
an obstacle, as we shall attempt to discuss in this work. We shall put emphasis 
on astrophysical tests of some versions of string theory entailing Lorentz-Invariance 
Violating space-time foam structures for the quantum gravity ground states. 

This brings us to the meaning of the word MAGIC used in this review. Here 
MAGIC is an acronym pertaining to the initials describing the full name of a Physics 
Instrument (M. A.G.I. C = .Major Atmospheric Gamma-ray Jmaging Cherenkov 
telescope). Specifically, it refers to a Telescope based on the Canary Islands obser- 
vatory (c.f. fig. 2), dedicated to the study of Cherenkov radiation emitted by highly 
energetic cosmic particles as they enter our atmosphere. From the study of the 
emitted Cherenkov radiation one can deduce several important conclusions on the 
nature of the initial particle and through this to try to understand the mechanisms 
of production of such energetic cosmic particles. 

On July 9th 2005, the telescope observed 3 (c.f. fig. 2) very high energy gamma 
rays from the active galactic nucleus Markarian 501 (Mkn 501), which lies at red- 
shift z — 0.034 (i.e. about half a million light years away from Earth). The highest 
energy photons with energy of order 10 TeV (1 TcV = 10 3 GeV = 10 12 cV), were 
delayed up to four minutes as compared with their lower-energy counterparts (in 
the 0.6 TeV or lower range) (c.f. fig. 3). This was the first observation of such a 
distinct delay. 

Three years after the MAGIC observations, in September 2008, the FERMI 
(formerly known as GLAST) Satellite Telescope 4 , also observed time delays of the 
higher-energy photons, from the distant Gamma Ray Burst (GRB) 080916c 5 , at 
red-shifts z = 4.35, and later on from GRB 090510 6 , at red-shift z = 0.9 and 
from GRB 09092B, at redhifts z = 1.822 7 . GRBs are cosmic explosions of titanic 
proportions, due to collapsing massive stars at distant parts of the universe. As we 
shall discuss in this review article, the pattern of these delays fits 8 a string model 
of quantum-gravity-induced refractive index, with the pertinent quantum gravity 
energy scale being essentially the same as that inferred from the MAGIC observa- 
tions (of order 10 18 GeV). Viewed as a lower bound, this scale is also compatible 
with that obtained from other Gamma-Ray data of the H.E.S.S. Collaboration 9,1 °. 
However, as we shall discuss in detail in this article, in order to accommodate the 
findings of MAGIC with the recently observed time delays of photons from the ex- 
tremely short Gamma-Ray burst GRB 090510 6 , provided that the latter can be 
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Fig. 1. The diagram shows the current (quite limited) understanding of the M-theory (interior 
of polygon graphs) unifying all known string theories, which appear as its low-energy limits. The 
left diagram shows the understanding we had in the mid 90s and is placed here for comparison 
of the progress made. The characterisation useless in the 1995 situation (upper panel) indicates 
the ability of the corresponding string theory to include the standard model group. The red- 
shaded interior regions lead to a partial understanding of M-theory issues, either through F-theory 
(termed F-th) considerations or rational conformal field theory (RCFT) analysis. The discovery 
(in 1996) of solitonic states in string theory, termed D(irichlet) brancs changed the picture. Many 
of the previously thought as useless string theories can incorporate standard model groups in 
non-trivial ways, involving intersecting brane situations (termed Inters, on the lower panel) and 
compactification to magnetised manifolds (termed Magn. in the right picture), orbifold (termed 
Orb.) or Calabi-Yau spaces (CY) etc. The M low-energy theory on the lower right corner of the 
super-unification graphs is identified now with the 11-dimensional supergravity (Pictures taken 
from Ref. 2 .) 

trusted a , one needs inhomogeneous densities of defects in the foam 11 , and in this 
sense such astrophysical observation are quite essential in falsifying models. 



a Indeed, there are uncertainties in this measurement concerning, for instance, the precise emission 
time of photons due to the duration of pre-cursors to burst activities and other such issues, which 
need to be confirmed by other measurements of similarly short bursts, that presently are lacking. 
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Fig. 2. Major Atmospheric Gamma-ray Imaging Cherenkov Telescope at the Canary Is- 
lands (Spain) Observatory (right panel). The telescope observed very high energy gamma 
rays, with energies up to the order of 10 TeV, from the active galactic nucleus Markar- 
ian 501 (radio image on left panel, by J.M. Wrobel and J.E. Konway, picture taken from 
http : / /www . vlba . nrao . edu/what is/mark . html). 

Of course, it goes without saying that the delay effects may be due to the con- 
ventional astrophysics of the active galactic nucleus or Gamma- Ray Burst (source 
effect), which, however, as we shall discuss below, is not well understood at present. 
In fact, currently there seem to be no consensus among the astrophysicists on the ap- 
propriate mechanism for the production of such high-energy photons at the source. 
These uncertainties prompted more ambitious, although admittedly far-fetched, ex- 
planations 12 , pertaining to new fundamental physics, affecting the photon propaga- 
tion due to space-time foamy vacuum structures that lead to modified dispersion re- 
lations for photons 13 - 15 ' 14 . If true, this would be a clear departure from the Lorentz 
invariant energy (E)-momentum (p) relations of Special Relativity, E — \p\c. 

The reader should bear in mind that at small length scales, of the order of the 

Planck length, lp = \J^%p- = jj^. ~ 1 0~ 35 m, which is the characteristic scale at 
which quantum gravity effects are expected to become dominant, the structure of 
space time may be quite different from what we perceive at our (low-energy) scales. 
It might be even discrete and non-commutative, that is the space-time coordinates 
(as we perceive them at present) might be average values of non commutative quan- 
tum operators. Moreover, one may have highly curved non-trivial fluctuations of the 
space-time metric, giving space time a "foamy" structure 16 . In such complicated 
Quantum Gravity (QG) vacua, the concept of Lorentz symmetry might break down 
at these short length scales, pointing towards the possibility of spontaneous Lorentz 
symmetry breaking by the QG ground state, since if Lorentz symmetry is intact, 
it strictly prohibits such modifications in the dispersion relations. Such departures 
from the standard Special Relativity form of dispersion relations correspond to the 
generation of a momentum dependent mass gap for the photon, and hence a non 
trivial refractive index, as the photon propagates through the medium of quantum 
gravity. 

It is the point of this review to touch upon such issues, through the description 
of the most important physical consequences of such a breaking. As we shall see, 
surprisingly enough, many models of QG that entail such a breaking can already 
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Fig. 3. The observations of the MAGIC telescope 3 regarding very high energy Gamma Rays 
(with energies in the TeV range) showed that the most energetic photons were delayed up to four 
minutes as compared with their lower-energy counterparts (in the 0.6 TeV or lower range). The 
figures show light curves (LC), i.e. the photon flux vs. time of arrival. Observations like this may 
be used to prompt new fundamental physics on the structure of space-time. The lower panel shows 
LC at different energy ranges, demonstrating clearly the time delay (of order of 4 ±1 minutes) of 
the more energetic photons. 



6 N.E. Mavromatos 



be falsified in current astro-physical experiments, which set very stringent bounds 
on Lorcntz symmetry breaking. Some of the models we are going to discuss have 
experimental consequences that can be tested (or are already falsified(?)) in Nature, 
at least in principle. In this article we shall discuss in some detail only one class 
of theories of QG that can entail such a breaking, which is a subset of the modern 
version of string theory, including brane (domain-wall-like) defects in space time. 
Such defects will play the role of the non trivial space time structures that would 
be deemed responsible for the spontaneous breaking of the Lorentz symmetry by 
the ground state of these systems. 

As we shall discuss here, there are very stringent conditions for such exotic 
explanations of the MAGIC observations to come into play in agreement with the 
plethora of many other astrophysical tests of Lorentz symmetry, existing currently. 
Specifically, the space-time foam must be transparent to electrically charged probes, 
such as electrons, while photons should exhibit non-trivial refractive indices in this 
theory, but with no birefringence effects. Moreover, the local effective lagrangian 
description of the effects, that is a low-energy representation of the QG medium 
effects in terms of local higher derivative operators in fiat space-time backgrounds, 
should break down. 

Thus, although one cannot exclude the possibility that both effects, source and 
propagation due to quantum gravity, may be simultaneously responsible for the 
observed photon delays in the MAGIC experiment, nevertheless the available the- 
oretical models that do the job are very limited. The thesis of this article will be 
that only certain models 17 > 18 > 19 of (the modern version of) string theory, includ- 
ing space-time defects, whose dynamics breaks Lorentz symmetry and provides a 
"foamy" structure of space-time, can satisfy the above conditions and thus offer 
an explanation for the MAGIC photon-arrival-times anomaly, consistently with all 
other current astrophysical constraints on Lorentz symmetry violations. The string 
modcl(s) provide an explanation for the observed photon time delays, in a natu- 
ral range of the string coupling and mass scale, and avoid all the other stringent 
constraints of Lorentz Invariance Violation coming from non observation of bire- 
fringence effects or of very high energy photons 20 . 

Our point in this article is not to advocate string theory as a superior candidate 
to other quantum gravity models available to date, but rather to discuss situations 
in which string theory predictions can be falsified by experiment. And high energy 
photon astrophysics may be a useful arena for this purpose! Of course, it goes 
without saying that, at present, we are very far from reaching any conclusions on 
such matters. Definite falsification of these stringy quantum gravity scenarios, if at 
all possible, would require a plethora of further studies by means of other high- 
energy astrophysics or particle physics processes and observations. 

The structure of the article is as follows: in the next section, 2, as a means 
of introduction, I discuss issues pertaining to the violation of Lorentz invariance 
in media or quantum electrodynamics vacua with non-trivial vacuum refractive 
index, e.g. the case of thermal plasma. Such cases may be thought of as (simplified) 
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analogues of Lorentz- Violating Quantum Gravity (QG) space-time foam vacua, but 
they also characterize the source regions of the cosmic high energy Gamma Rays. 
In this latter sense, it is important to understand photon propagation in those cases 
first, so as to disentangle possible source effects from the QG-induced ones. In section 
3, 1 proceed to the main topic of this review, namely a discussion on the string/brany 
model of space-time foam. The model involves membrane-like defects in a higher- 
than-four dimensional space time, with our world being viewed as a hyper-membrane 
(D(irichlet)-brane : ) embedded in this space time. The model is of the kind of large- 
extra-dimension models to be tested at LHC and future colliders, with the important 
ingredient of having point-like space-time (D-particle) defects, responsible for the 
"foamy" structure of space time. I describe first the cosmological features of this 
model, specifically in connection with modifications of the energy budget of the 
Dark Sector (Dark Matter and Dark energy), compared to the standard ACDM 
model. Then, I proceed to discuss its "optical properties", namely the induction 
of effective vacuum refractive indices, the absence of birefringence, and the break 
down of a local effective lagrangian formalism insofar as the description of the 
induced time delays of high energy photons (due to interactions with the foam) are 
concerned. I also discuss a novel ultra-violet "cutoff" that may characterize ultra 
high energy photon scattering for specifically stringy reasons. In section 4, I discuss 
high energy Gamma Ray astrophysics as a potential probe of the stringy space-time 
foam model. First, I analyze possible interpretations of the MAGIC (and FERMI) 
observations, including exotic ones involving QG dispersive media. Then, I discuss 
bounds and sensitivities of various astrophysical experiments, and state carefully 
the stringent requirements that must be met by a theoretical model of QG, in order 
for the observed delays in the MAGIC and FERMI Telescopes to be attributed to 
effects due to a space-time-foam medium, in agreement with all other current tests. 
All these requirements are surprisingly respected by the string model discussed in 
section 3. Finally, section 5 contains concluding remarks and an outlook. 

A note is in order at this point concerning the units used in this work: throughout 
the article, unless otherwise stated, we shall work in natural Planck units, in which 
h = c = 1. In these units, length and time are identified, and they are inversely 
proportional to mass or energy. The latter are also identified and expressed in units 
of multiples of eV (=1.6 x 1CT 19 Joules), in particular GeV (=10 9 cV) and TcV 
(=10 12 eV) in this work. From time to time, for concreteness, the speed of light 
in vacuo c and Planck's constant h or h = h/2n may appear explicitly in some 
formulae. 

2. Lorentz Invariance and the Vacuum Structure of Quantum 
Fields 

One of the cornerstones of Modern Physics is Einstein's theory of Special Relativity 
(SR), which is based on the assumption that the speed of light in vacuo c is an 
invariant under all observers. In fact, this implies invariance of the physical laws 
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under the Lorentz transformations in fiat space times, and the role of c as a universal 
limiting velocity for all particle species. 

The generalization (by Einstein) of SR to include curved space times, that is the 
theory of General Relativity (GR), encompasses SR locally in the sense of the strong 
form of the equivalence principle. According to it, at every space-time point, in an 
arbitrary gravitational field, it is possible to choose a locally inertial ('free-float') 
coordinate frame, such that within a sufficiently small region of space and time 
around the point in question, the laws of Nature are described by special relativity, 
i.e. are of the same form as in unaccelerated Cartesian coordinate frames in the 
absence of Gravitation. In other words, locally one can always make a coordinate 
transformation such that the space time looks flat. This is not true globally, of 
course, and this is why GR is a more general theory to describe gravitation. The 
equivalence principle relies on another fundamental invariance of GR, that of general 
coordinate, that is the invariance of the gravitational action under arbitrary changes 
of coordinates. This allows GR to be expressed in a generally covariant form. 

In such a locally Lorentz-invariant vacuum, the photon dispersion relation, that 
is a local in space-time relation between the photon's four-wavevector components 
k^ = (w, k) (where to denotes the frequency, and k the momentum) reads in a 
covariant notation: 

Wri^ = (1) 

where repeated indices /x, v = 0, 1, ... 3, with referring to temporal components, 
denote summation and n^ denotes the Minkowski space-time metric, with compo- 
nents 7/00 = — lj Voi — Vio — 0, r/ij = r]ji = Sij , i = 1, 2, 3 with Sij the Kronecker 
delta symbol. 

The above relation (1) implies the equality of all three kinds of photon velocities 
in vacuo that stem from its wave nature (due to the particle- wave duality relation): 

phase : v ph = ^ = ^ = c 

group : v gr = Jg- = = c , n gr (w) = n(w) + uj^T 

front : f front = c/n(oo) = c (2) 

since the phase and group refractive indices of the trivial vacuum equal unity 
n(w) = n gr (uj) = 1. For brevity we shall work from now on in units where c = 1. 

2.1. Photon Propagation in Conventional Media and in 
Non-trivial field-theory vacua with a refractive index 

The above results (2) change significantly when light propagates in a material 
medium, in which its speed is different from c in vacuo. This is due to the non 
trivial refractive index the material has, as a result of the electromagnetic interac- 
tions of the photon with the electrons in the medium. As we shall discuss later on, 
this case seems to bear some quite instructive analogies with our string model of 
space time foam, which lies at the focus of our attention in this review. 



String Quantum Gravity, LIV and Gamma-Ray Astronomy 9 



m 



The simple model of quantum oscillators has been adopted by Feynman 21 as 
a simplified but well motivated analogue for describing the situation in the case of 
photons in ordinary media. In that case, the electrons of the medium, of mass m, are 
represented by simple harmonic oscillators, with frequency cj , which provides the 
necessary restoring force during the scattering of light off the electrons in order to 
keep the latter oscillating around their initial position. In that problem the induced 
refractive index is obtained by the reduction of the phase velocity of the photon 
wave. 

We consider 21 the electrons, of mass m, as forced simple-harmonic oscillators 
with a resonant frequency uj , subject to the force F exerted by an oscillating 
external electric field of frequency u>: F = eEoe tult , where e is the electron charge. 
The corresponding equation of motion is: 

(JL x + bJ 2 x ^ =eEo e i " t (3) 

Assuming for concreteness a plate of thickness Az, representing the medium through 
which an electromagnetic wave travels, where z is perpendicular to x, one can 
compute in a standard way the electric field E a produced by the excited atoms: 

F a = -^i , e 2° ^ e*^*-*), (4) 

where e is the dielectric constant of the vacuum and n e is the area density of 
electrons in the medium (plate), which is given by n e = p e Az, where p e is the 
volume density of electrons. 

We next recall that light propagates through a medium with a refractive index 
n with a speed c/n, causing a delay At in traversing the distance Az, given by: 

At={n-l)Az/c . (5) 

Representing the electric field before and after passing through the plate as 
^before = £ e J " (t - z/c) and £ aftcr = E e LUl{ - t - z / c - {n - 1 ^ z l c \ in the case of small de- 
viations from the vacuum refractive index we have: -E a ftcr — EQe lLU ^ t ' z ^ c ^ — i[uj(n — 
l)A2/c]£'oe l " ( -* _z / c - ) . The last term on the right-hand-side of this relation is just the 
field E a produced in the region of space after the plate by the oscillating electrons. 
We then obtain from (4): 

(»- 1 )^=554fc»J' (6) 
and hence the following formula for the refractive index in a conventional medium: 

n = 1 + 9 p '"f 2V (7) 

We see in (7) that the refractive index in an ordinary medium is inversely propor- 
tional to (the square of) the frequency to of light, as long as it smaller than the 
oscillator frequency, where the refractive index diverges. 

If the couplings of the two polarizations of the photon to the electrons in the 
medium are different, the phenomenon of birefringence emerges, namely different 
refractive indices for the two polarizations. Moreover, we see from (7) that the 
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propagation of light is subluminal if the frequency (energy) of the photon to < uj , 
whereas it is super luminal for higher frequencies (energies) 21 . This reflects the fact 
that the phase shift induced for the scattered light can be either positive or negative, 
but there such a superluminal refractive index causes no issue with causality, since 
the speed at which information may be sent is still subluminal. 

As we shall see later on, this conventional model will provide us with a rather 
good analogue of what happens in some string models of quantum space time foam, 
also characterized by non-trivial refractive indices, which we shall analyse in section 
3. However, as we shall see there, contrary to the conventional situation discussed 
in this section, in the string case the refractive index is found proportional to the 
photon frequency, while the effective mass scale that suppresses the effect is the 
quantum gravity (string) scale and not the electron mass as in (7). 

A word of caution on this point concerns the disentanglement of the time delay 
due to a non-trivial refractive index (5) from a modified dispersion relation. In some 
cases the two are equivalent. In the string model, however, discussed in section 3.4, 
the pertinent delays will be associated with stringy uncertainties of the foam, and 
will not be directly related to modified dispersion relations for photons 18 > 22 . 

The role of a non-trivial refractive index material can be played under certain cir- 
cumstances by a non-trivial vacuum in which photons propagate, such as quantum 
electrodynamics at finite temperature plasmas 23 or the Casimir vacuum between 
parallel capacitor plates 24 (or other geometries, as long as the space is bounded 
appropriately) (c./. fig. 4). In such cases, the loop corrections due to vacuum polar- 
ization, i.e. creation and annihilation of virtual electron-positron pairs, in quantum 
electrodynamics (QED) result in a modified photon propagator, and a non-trivial 
group velocity and refractive index n(w) ^ 1. The reason for this is the breaking 
of Lorcntz invariance, due to either the existence of spatial boundaries (Casimir 
vacuum) or finite temperature (thermal vacuum in case of plasmas). Moreover, one 
may consider quantum electrodynamics in a homogeneous and isotropic Friedman- 
Robertson- Walker (FRW) expanding- Universe background 25 and examine the non- 
trivial effects of vacuum polarization on photon propagation there. This was in fact 
historically the first instance where the effects of curvature induced superluminal 
propagation for low-frequency photon modes. 

The situation concerning all the above cases can in fact be represented in a 
rather unified way by the following formula for the group velocity of (low-energy) 
photons, due to vacuum polarization in the context of four-dimensional QED in 
such non- trivial vacua 23 : 

V = 1 - ^ 2 lk ■ ^ 
This formula is valid in all cases except the Gravitational Background case of 25 , 
where a 2 should be replaced by am 2 GN, with Gn the Newton gravitational con- 
stant. Here, p denotes the energy density relative to the standard vacuum, and can 
be negative or positive. In the (four space-time dimensional) Casimir vacuum, for 
instance, the presence of boundaries in space (e.g. capacitor parallel plates at a 
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Fig. 4. The Casimir Vacuum of Quantum Electrodynamics: in the compact space between the 
parallel plates (or more generally in other compact geometries), the quantum fluctuations of the 
electromagnetic field are responsible for the macroscopically measured Casimir force between the 
plates, scaling with the size L of the compact space as L~ 4 in four space-time dimensions. The 
force can be either attractive or repulsive, depending on the geometric set up. However, besides this 
force from "nothing" , the Casimir Vacuum provides another interesting phenomenon in quantum 
physics. According to calculations by Scharnhorst and collaborators, the presence of the boundaries 
breaks manifestly Lorentz invariance, and results in modified dispersion relation, and hence a non- 
trivial refractive index, for the virtual photons of the non-trivial Quantum Electrodynamics vacuum 
in the compact region. From this latter perspective, the situation is entirely analogous to what 
happens in thermal plasmas, such as the ones occurring in the interior of stars or source regions 
of active galactic nuclei, which we are interested in in this work. There, the periodic boundary 
conditions of the Casimir case are replaced by similar ones due to the finite temperature, which 
also break Lorentz symmetry of the ground state. 



distance L in the simplest geometry), imply loss of low-energy photon modes, and 
as such the energy density of the vacuum is lowered relative to the standard one 
(without the boundaries). In this case p = — 72 q L4 < in (8), and one has the 
phenomenon that the effective low-energy photon modes appear to propagate in a 
superluminal way, v gr > 1 (we give here the formula for propagation perpendicular 
to the plates for concreteness and ease of comparison with the plasma case later 
on): 



•v. 



Casimir _ i , ll7r 2 „,2 1 i (n\ 



As explained in 23 , there is no contradiction with relativity here, as this result applies 
only to low energy photons, with energies much lower than m e , and indicates the 
effective loss of degrees of freedom due to the spatial boundaries. 

In the plasma case at finite temperature T, which characterizes the interior 
regions of stars or active galactic nuclei, of interest to us in this work, there is a 
formal analogy 23 between the roles played by the temperature T and the plate 
separation L in the Casimir vacuum set up. They both break Lorentz invariance in 
a rather similar way. In fact, there is a correspondence between the Casimir and 
thermal-plasma vacuum formulae by replacing IT by L . In the plasma case, the 
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low energy photon modes with momentum k 2 <C m\ (k = \k\) have a group velocity: 
VgI (kT« m t) = l + %£a 2 (%.y > 1 (10) 
which can be directly compared with the Casimir one (9) upon the substitution 

For future reference we give the group velocity for the high-energy (as compared 
to the effective QED scale m c ) photon modes 23 in the finite T plasma case: 

V(^>> m e 2 ) = < 1 (11) 

which again can be compared with the corresponding Casimir vacuum case (9), 
for photon modes perpendicular to the plates 24 , upon the replacement 2t — > 
The reader should notice here the momentum- (A:) dependence of the subluminal 
velocity for the high-energy photon modes. In such non-trivial vacua, therefore, high 
energy photon modes will exhibit a momentum dependent subluminal refractive 
index, which will affect their arrival times at the observation point, if one considers 
simultaneous emission of modes within a certain energy range. From (11) we observe 
that the higher the momentum k of the high-energy photon mode the higher the 
group velocity, since 9, v( fc ^ >rn c) _ ^hi(i^)^ > o, and thus fast modes will 

arrive first if emitted simultaneously in this vacuum. Similar effects also characterise 
the curved background case of 25 , where again the high- momentum photon modes 
are subluminal. 

We remark at this point that the relation (11) can characterise the source regions 
of active galactic nuclei, and thus such effects can be responsible for inducing time 
delays of photons from these regions. We shall come back to this point when we 
shall discuss the MAGIC delays in section 4.3 below. As we shall see there, however, 
the plasma-induced time delays are much smaller than the observed delays in the 
MAGIC experiment, and hence cannot provide the dominant explanation for the 
effect. 

At this point it is instructive to make some clarifications regarding the speed 
of light and causality, that is the fact that signals do not arrive before they occur. 
The phase, front and group velocities have all been found to exceed the value of the 
speed of light in vacuo. However, this is not in conflict with causality. 

Indeed, the various light velocities do not have to be subluminal, as they carry 
no information. Information can be transmitted by (or not) sending pulses. The 
information then appears at first sight to propagate with the group velocity, i.e. 
the velocity of the peaks of the pulses. However, as experimentally demonstrated 
in several instances since the 1980's, the group velocity of the photons can also 
be super-luminal. This is because the group velocity can be disentangled from the 
signal velocity, i.e. the speed by which information can be transferred, and thus 
it can exceed the speed of light in vacuo, c, without contradicting the causality 
requirements of Special Relativity. Group velocities larger than c can occur e.g. in 
tunneling experiments and appear to lead to superluminal transmission. However, 
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there is no contradiction with causality or Special Relativity in such cases. The error 
lies in identifying the peak of the pulse with the temporal position of the carried 
information. For example, a Gaussian-shaped pulse can be detected long before its 
peak due to the rise of intensity at earlier times. Therefore, a different kind of signal 
must be considered, where no information at all is sent out before a certain moment 
of time. For such signals, it can be proven that the earliest time at which that 
switching event can be observed is limited exactly by propagation with the vacuum 
velocity c. A so-called precursor is traveling with that speed, but is normally too 
weak to be detected, except in certain circumstances. For direct measurements of 
optical precursors in regions with anomalous dispersion the reader is referred, for 
instance, to ref. 26 . 

Super-luminal group photon velocities have been measured in laser pulses pass- 
ing through specially prepared materials in ref. 27 . Using special set up involving 
optical fibers, a super-luminal group velocity of photons (in the fiber) has been 
measured and found different from the signal velocity, which is defined as the speed 
of the front of a square pulse. It is the signal speed that cannot exceed c, due 
to causality, as mentioned above. The precursors (or "forerunners") of the signal, 
mentioned above, which travel with sub-luminal speed, arrive first before the main 
front. This signal velocity was measured for the first time directly in an experiment 
in 27 , and indeed was found to be less than c. 

This is an important issue that the reader should have in mind, especially when 
we discuss non-trivial vacuum refractive indices in several non trivial ground states 
of quantum systems, including gravity. 

In fact, as argued in ref. 28 (specifically for the case of Casimir vacuum 24 but 
the discussion can be generalised), super-luminal group velocities are "benign" as 
far as causality and compatibility with the kinematics of Special Relativity are con- 
cerned. In particular, the kinematics of Special Relativity requires only an invariant 
speed not actually a maximum one. Moreover, causality can be guaranteed in such 
super-luminal cases because the pertinent kinematics is equivalent to an "effective" 
metric in space-time, g^g, \i, v = 0, . . . 3, different from the Minkowski one, that de- 
scribes the kinematics in the non-trivial (Casimir) vacuum. The modified dispersion 
relations leading to super-luminal group velocities (9) acquire the form 



with the (inverse) Minkowski metric, a unit (space-like) vector orthogonal 
to the plates of fig. 4, and £ = 4050^ ° m 4 ' w ^ tn ^ tne distance between the Casimir 
plates. 

The basic point of the discussion in 28 , and how causality is maintained, is 
that the presence of the effective metric (13) in (12) widens slightly (due to the 
deviations of order £ <C 1 of the metric from the Minkowski one) the light cone 




1,2,3 



(12) 



(13) 
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in the direction orthogonal to the plates, and hence light in that direction travels 
at a speed ci; g ht > c, while light in the direction parallel to the plates still travels 
with speed c. In a given inertial reference frame, moving with four-velocity 
with respect to the rest frame of the apparatus of fig. 4, the photons inside the 
Casimir cavity travel, at a direction perpendicular to the plates, with a definite 
speed c[jgk t > c, which has only one value for each observer and is not universal 
among observers. In this way super-luminal group velocities (9) are compatible with 
causality, since violation of the latter occurs only if signals travel with the same 
speed greater than c in two different frames. 

As a final comment on properties of the speed of light in non-trivial vacua 
we mention some recent work in 29 according to which the above-mentioned fea- 
tures of superluminal phase-, and hence front-, velocity of light might have some 
effects on micro-causality, which is a property related to the local commutativity 
(or anti-commutativity in the case of fermions) of fields (i.e. the requirement that 
the vacuum expectation value of the commutator of two field operators at space 
like separations of x, y, (0\[A(x) , A(y)\0) be zero). This is a fundamental axiom of 
local quantum field theories. In fact, the authors of 29 argue that such correlators 
no longer vanish in the above-mentioned non-trivial vacua. 

Indeed, although the subluminal nature of the high frequency photon modes 
guarantees macro- causality, that is the property that signals are not measured before 
they arrive, nevertheless the superluminal nature of the phase velocity of the low- 
energy photon modes might 29 create some problems with the micro-causality in 
curved space time backgrounds. Microcausality is usually studied by examining 
the validity of Kramers-Kronig (KK) dispersion relation relating the (complex in 
general) phase refractive index at high frequency to that at low frequency, namely 

Ren(oo) - Ren(O) = - f / °° ^Imn( W ) (14) 

It should be remarked at this point that in standard optics this relation is assumed 
to be valid, but this may not be the case of (quantum) gravitational backgrounds, 
for which the situation is still far from being understood in detail, due to the lack 
of a fundamental underlying theory. 

The imaginary parts of the refractive index in (14) refer to absorptive prop- 
erties of the medium. In fact the KK dispersion relation follows from analytic- 
ity requirements of the (complex) refractive index in the upper half plane. If one 
assumes (14) for classical gravity, then it becomes evident that, by assuming a 
positive Imn(w) > 0, which is the standard property of an absorptive medium 
and follows from unitarity arguments of a quantum field theory, one arrives at: 
Ren(oo) = l/u p h(oo) > Ren(0) = l/u p h(0), i.e. v p h(oo) > w p h(0). Hence, a superlu- 
minal low-frequency phase velocity, as seems to be the result of vacuum polarization 
in QED in curved space times 25 , would imply a superluminal front velocity, with the 
consequent violation of causality. Assuming the KK dispersion relation, therefore, 
the only way out of this would be to assume a negative Imn(w) < for gravity, i.e. 
classical gravity does not behave as an absorptive but rather as a gaining medium. 
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In nature one can meet situations like this, for instance in laser-atom interactions, 
which can induce gain, leading to superluminal low frequency phase velocities, but 
preserving a front velocity equal to the speed of light in vacuo c and hence the KK 
dispersion relation. 

However, as indicated in 29 , QED seems to violate the KK relation itself, by in- 
troducing non analyticities in the corresponding refractive index. Since this relation 
relies only on analyticity of the refractive index in the upper half plane, a possible 
violation would indicate the breakdown of micro-causality. 

However, the alert reader should bear in mind that the requirement of micro- 
causality actually requires Lorentz invariance of the Scattering matrix, and hence 
it may not be valid in a generic curved space time theory, where Lorentz invariance 
may be violated. In fact, as argued in 29 , it appears that QED violates analyticity in 
the presence of curved space times, and in this way KK relation fails, and causality 
is maintained. The authors of 29 found that in the high-frequency limit, the phase 
velocity always approaches c, so they determined Vf ront /c = 1. Moreover, they have 
also shown that where the background gravitational field induces pair creation, 
7 — > e+e~ the Imn(w) is indeed positive as required by unitarity. However, the 
refractive index n(oj) is not analytic in the upper half-plane, and the KK dispersion 
relation is modified accordingly. In particular, for the specific cases of space-time 
backgrounds of constant scalar curvatures R, they examined for concreteness, they 
found that there are branch cuts in the upper half plane of the refractive index, due 
to vacuum polarization (computed using the world-line formalism of QED in order 
to deal properly with the high frequency regime), which modify the KK relation 
such that, for conformally flat backgrounds, for instance, one obtains: 

Ren(co) - Rcn(O) = ^ > (15) 

One might think that this implies a violation of micro-causality. In our opinion, 
however, the issue of the refractive index of background gravity is not the end of 
the story in a dynamical space time, where back reaction effects of the propagating 
matter have to be properly taken into account, in particular at high energies. In this 
respect, there may be problems with the KK dispersion relation, associated with the 
fact that at the high frequency branch, uj — > oo, one is actually dealing with modes 
with super-Planckian energies, i.e. regimes where the full quantum gravity theory 
may come into play. Hence, it may be misleading to attempt to define concepts from 
flat space time at such regimes. In particular, back reaction effects of transplanckian 
modes may affect the unitarity properties of a quantum field theory. At such regimes, 
the structure of space time is highly curved, or even discrete, so any attempt to rely 
on analyticity properties of the scattering matrix may breakdown. 

Above the energy scales Mqq, where quantum gravity may set in, the struc- 
ture of the theory may change completely. Concepts like those of a local effective 
lagrangian or quantum coherence may be in jeopardy 30 ' 31 ; leading to violations of 
fundamental principles of ordinary quantum field theories, such as invariance of the 
effective lagrangian or correlation functions under CPT, i.e. the successive opera- 
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tions (at any order) of Charge conjugation (C), Parity (reflexion) and Time reversal 
(T) symmetries. Hence, there is no guarantee that the formal KK relation, which 
involves an integration all the way up to to — > oo, is valid above such high scales, 
in the complete situation, where back reaction effects are taken into account and 
gravity is no longer treated classically as a background. 

The issue of the refractive index of gravity depends crucially also on the type 
of the underlying microscopic theory. The reader should have these in mind when 
studying space-time foam theories of quantum gravity. In this review we shall dis- 
cuss in section 3.4 a particular type of string foam, leading to a subluminal refractive 
index, whose anomalous effects are proportional to the photon frequency. This is 
a result of non-trivial quantum interactions of photon with space-time stringy de- 
fects. In this model, the above issues on possible violation of micro-causality do not 
arise. Nevertheless, the effects of the stringy foam on low-energy particles cannot 
be described within the framework of local effective Lagrangians. 

2.2. Non-trivial Optical properties of the Quantum- Gravity 
Vacuum ? 

A truly unspeakable feature on the speed of light may appear when considering the 
ground state of Quantum Gravity per se as a non-trivial vacuum with non-standard 
optical properties, leading to a non-trivial refractive index. This may characterize, 
for instance, certain approaches in which path-integration over microscopic singular 
fluctuations of the metric field, such as Planck size black holes and other topolog- 
ically non-trivial configurations, implies a sort of "foamy" structure of the space 
time at small (Planck-size) length scales, over which photons can propagate. This 
proposal was made initially by J. A. Wheeler 16 . In such situations, the concept of 
a local Effective Quantum- Field Theory Lagrangian may break down, and the situ- 
ation resembles that of a quantum-decoherent motion of matter in open quantum 
mechanical systems interacting with an environment 30 > 31 > 32 . 

The ground state of such QG foam situations may behave as a (subluminal) 
refractive medium, as suggested originally in 3 , which to our knowledge, constitutes 
the first concrete attempt to consider the non-trivial optics effects effect of such 
vacua on massless particle (photon) propagation, based on earlier works on non- 
critical strings in black hole backgrounds 34 . We note at this stage that similar 
ideas on modified dispersion relations for particles in a quantum gravity medium, 
but on a purely phenomenological basis, without any attempt to present concrete 
models, have also been advocated in ref. 35 . Subsequent to the suggestion of 33 , 
modified dispersion relations have also been suggested to characterize some Loop 
quantum gravity ground states 36 . 

In the quantum gravity case, the effective modified dispersion relations of the 
low-energy matter theory, assumes the generic form: 

oo 

E> = \pf + m 2 + 5>|pf (JJL)" (16) 

71=1 
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where m is the rest mass of the probe, and c„ are constant coefficients, with signature 
and values that depend crucially on the type of theory considered. It is not clear 
whether the series converges, or is resummable, as this information depends crucially 
on the details of the underlying microscopic theory. The important point to notice 
in (16) is that the natural suppression scale of the (Lorentz-symmetry-violating) 
correction terms of the standard Special Relativity dispersion relation is that of 
Quantum Gravity itself, Mqq. According to what was mentioned earlier, there is 
no fundamental issue with violation of (micro-)causality 28 or contradiction with 
the principles of Special Relativity, at least for some cases where super-luminal 
group velocities arise from QG anomalous dispersion relations (16). The important 
point to observe is that, as in the Casimir vacuum case (12), effective metrics can 
be found, at least in some cases, which describe the QG anomalous dispersion 
(16). Thus, causality can be saved for those QG cases by applying the same logic 
as for the Casimir vacuum, discussed previously. Nevertheless, as we shall discuss 
later on in subsection 4.6, in cases where super-luminal modes exist, one has QG- 
induced birefringence phenomena, which are severely constrained by the current 
phenomenology (i.e by the absence of the relevant signals). 

Before closing this subsection we also mention that there is another approach 
towards modified dispersion relations of particles, of the type (16), based on the 
so-called doubly special or deformed special relativity (DSR) theories 37 > 38 . In such 
approaches, which are formulated on fiat space times, the modification to the dis- 
persion relation arises by the postulate that the local symmetry group of space 
time is no longer the Lorentz group but a different one. In their original version 37 , 
DSR modified dispersion relations were obtained by the requirement that the length 
scale of quantum gravity ("Planck") remain invariant under transformations, which 
thus leads to deviations from the Lorentz group. There is no unique prescription 
to achieve this, however, and in this way one may even arrive at models where 
there is no upper limit in momentum. In their subsequent version 38 , DSR models 
postulated the existence of upper limits in velocities of species, and in this way the 
modified local group was determined by appropriate combinations of dilatations 
and translations. However not all DSR theories exhibit modified dispersion 39 , and 
the issue there lies on what one defines as the physical momenta. 

It is unclear to us whether such theories are fundamental or effective, and in 
the former case, whether they can be quantized consistently. For instance, a known 
problem is the behaviour of multi-particle states in such DSR fundamental theories. 
Another problem for quantizing DSR theories with modified dispersion relations is 
associated with their lack of locality, as pointed out recently in , although currently 
there is an ongoing debate as to whether there is a resolution to this problem. 

We shall not discuss these theories further here, apart from mentioning that, 
as in the stringy model we shall consider in section 3 33 > 41 > 4 2,43,i8,i9^ they are not 
characterised by birefringence effects 44 , leading only to subluminal propagation 
(defined appropriately 37 ' 38 ). There is, however, an important difference from the 
stringy model: the action of gravity is universal among particle species in DSR 
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theories, and therefore the resulting modification in the dispersion relation pertain 
to all species. In contrast, for specifically stringy reasons to be outlined in section 
3.4, in the string model only photons and at most electrically neutral particles 45 are 
allowed to interact non-trivially with the stringy defects, and are thus subjected to 
non-trivial modification of their dispersion relations. Moreover, in our string foam 
models, the modified dispersion relations arise as a result of local distortions of 
space time due to the recoil movement of space time defects, inducing metrics that 
depend on momenta (Finsler-like 46 ). The underlying string model is perfectly well 
defined as a quantum theory, and the associated effects are perfectly consistent 
with stringy uncertainty principles, as we shall discuss below. Hence, criticisms of 
the type appearing in 40 for DSR are not relevant for our string foam models. 
Nevertheless, as we shall also discuss here, there seems to be a breakdown of a local 
effective field theory description of the effects of string foam on the propagation of 
particles, which in fact proves crucial for the model to be able to evade 22 stringent 
bounds on Lorentz Violation from ultra high energy cosmic ray physics 20 . 

3. Lorentz-Invariance Violating Stringy Space-time "Foam" 

In string theory, the simplest way to obtain Lorentz-Symmetry Violation, consis- 
tently with the conformal properties of a c-model vacuum, as appropriate for a 
first-quantized version of strings, is to consider constant electric and/or magnetic 
background fields. The Lorentz symmetry violation is associated with the preferred 
reference frame specified by the direction of the background field. In a er-model lan- 
guage, the presence of a constant electric or magnetic background field, in which a 
string propagates, is described by deforming the world-sheet action by appropriate 
antisymmetric tensor = —B vil backgrounds. Constant Magnetic fields are asso- 
ciated with the spatial components of the B-field, B^ ^ 0, while a constant electric 
field is associated with the Oi-component B^ ^ 0, with the temporal component 
and i the spatial component in the direction of the electric field. 

The presence of a B-ficld background leads also to non commutativity (NC) of the 
target-space string coordinates 47 . For spatial- coordinate NC one needs only spatial 
components of the £>-field to be non zero (corresponding to constant "magnetic" 
fields in the simplest case), while for time-space non commutativity one needs 48 
B M 7^ 0, which is equivalent to constant "electric" background fields, in its simplest 
version. 

3.1. Strings in Constant Background Fields: the simplest case of 
Lorentz Violation 

To understand better the above points, we should note that, in general, the _B-field 
tr-model vertex operator on the world sheet X in which open strings propagate has 



the form 




(17) 
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where B^ v = — a,(3 = 1,2 arc world-sheet indices, (i,v = 0,1, . . . D arc target 
space indices and e a/3 is the two-dimensional antisymmetric (Levi-Civita) symbol. 
In the modern version of strings, including branes 1 , the dimensionality D may 
represent the target-space dimensionality of branc worlds, which the ends of the 
open strings are attached to. In the case of a constant background gauge field, 
with constant field strength F$) one identifies B^ v = —^F$), as follows by the 
application of Stokes theorem to the vertex operator of the corresponding gauge 
potential 

/ A^X^ = [ e a ?d a (A^dpX^) = -i J E e^F^dvX^dpX" , 

where 9S the boundary of the world-sheet and = d^A v — d^A^ . 

The presence of the B-field leads to mixed-type boundary conditions for open 
strings on the boundary &D of world-sheet surfaces with the topology of a disc, we 
consider here for concreteness and relevance to our discussion below: 

9„ v d n X v + B^drX^gv = , (18) 

where <7 M „ denotes the metric of target space-time. 

In what follows, we shall concentrate for concreteness to the case of constant 
electric background fields E. This will be relevant for our subsequent discussion on 
space-time foam. Considering commutation relations among the coordinates of the 
first quantized er-model in the above background, one obtains a non-commutative 
(NC) time-space relation 48 , with the pertinent spatial coordinates along the di- 
rection of the electric field, say X 1 for concreteness and brevity. The NC relations 
assume the form: 

[X\t]=i6 w , 6 ai (=-9 w )^9 = 1 L T M^ (19) 

where t is the target time. The quantity Ei = jjf- and E c = -^p- is the Born-Infeld 
critical field. The space-time NC relations (19) are consistent with the space-time 
string uncertainty principle 49 

AX At > a' (20) 

As discussed in detail in refs. 48 > 47 ; there is also an induced open-string effective 
target- space-time metric. To find it, one should consider the world-sheet propaga- 
tor on the disc {X^{z,~z)X u (0, 0)), with the boundary conditions (18). Upon us- 
ing a conformal mapping of the disc onto the upper half plane with the real axis 
(parametrized by r £ R) as its boundary 47 , one then obtains: 

<X"(r)X"(0)) = -a'Cen, dcctnclnr 2 + i^e(r) (21) 

where e(r) is the function that is 1 (—1) for positive (negative) r. The non- 
commutativity parameters Q^ v are given by (19), while the effective open-string 
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metric, due to the presence of the electric field E, whose direction breaks target- 
space Lorentz invariance, is given by: 

ff °P en ' electric = (l - Ef) v, , u,v = 0,l 

ff opc„,clcctric = r]iiv il v = aU Qther valueg ) (22) 

Moreover, there is a modified effective string coupling 47 < 48 : 

1/2 



- ,9s (l - £ 2 ) (23) 



Notice that the presence of the critical "electric" field is associated with a singu- 
larity of both the effective metric and the non commutativity parameter, while the 
effective string coupling vanishes in that limit. This reflects the destabilization of 
the vacuum when the "electric" field intensity approaches the critical value, which 
was noted in 50 . 



3.2. Non-commutative Effective Field Theory 

Once non-commutativity is established, one may proceed to write down a generic 
(low-energy) effective field theory, by following the rather generic approach of 51 . For 
constant electric fields, which we are restricting our attention to for the purposes of 
this review, the resulting effective field theory is of a type that can be accommodated 
within the so-called Standard Model Extension (SME) framework 52 . 

For brevity, and relevance to our phenomcnological discussion in this work, we 
shall restrict ourselves here to a discussion of Non-commutative Quantum Electrody- 
namics (NCQED), following 51 , for a generic non-commutativity among coordinates 
of the form 

[x», x v \ = ie» v , e^e^ > o , (24) 

which is certainly of the type induced by a constant electric field (c.f. (19)). The 
condition 9^9^ > guarantees perturbative unitarity of the effective theory 51 . 

One approach to constructing a non-commutative (NC) quantum field theory is 
to promote an established ordinary theory to a noncommutative one by replacing 
ordinary fields with noncommutative fields and ordinary products with Moyal ★ 
products, defined by 

f*g{x) = cM¥^d x »d v »).f(x)g(y)\ x=y . (25) 

For gauge theories, one can define appropriate extensions of the pertinent gauge 
tansformations, and the covariant derivatives acting on the NC fields. Denoting by 
a caret symbol over a field the corresponding NC extension, we write for the gauge 
potential, A^, the fermion field ip, the Maxwell tensor — d^A^—d^A^ — ilA^, A v \ 
and the corresponding NC extension of the covariant derivative: D^ip = d^ip — iA^* 
'ip. The gauge invariant NCQED Lagrangian then, proposed in 51 , reads: 
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with /* g = f * D^g — D^f -kg. Any NC field theory, and hence (26), violates 
Lorentz symmetry, since the NC parameter 9^ carries space-time indices. In the 
case of strings in external background hclds the violation of Lorentz symmetry by 
the direction of the field is obvious. 

An important point is in order here, concerning the charge assignments of (26). 
From the definition of the NC gauge transformations and the associated gauge 
covariant NC derivatives only particles with charges q = 0,±1 are allowed, which 
certainly seems problematic from the point of view of associating (26) to an effective 
action of string theory in constant external electric field backgrounds, which in 
general incorporates the Standard Model. 

This issue is resolved 51 by the so-called Seiberg-Witten map 47 , which in the 
case of NCQED and for small NC parameter, we assume here, defines the physical 
excitations associated with the caret fields in (26) as follows: 



This leading-order form suffices for many purposes, since any physical non- 
commutativity in nature must be small. In fact, as we shall discuss later on in 
section 3.6, in our stringy space-time foam model the role of the non commutativity 
parameter is played by the recoil velocity of the heavy D-particle defects in the 
space-time foam, which is small. 

Plaguing the above expressions (27) to the action (26) and expanding in the 
small NC parameter 9^" leads to the following effective action that is physically 
equivalent to noncommutative QED to leading order in 9^ 51 : 



In this equation, one has redefined the gauge field — > qA^ to display the charge 
coupling of the physical fermion, and = d^ijj — iqA^ip as usual. This action is 
manifestly CPT invariant, and as such it constitutes only part of the SME. 

Let us now connect these results to those coming from strings in external elec- 
tric field backgrounds, which is our main topic of discussion here. In the context 
of strings in external background fields the above action can be obtained from the 
standard procedure of considering the equations of motions for the various back- 
ground fields, obtained from the requirement of world-sheet conformal invariance. 
Let us ignore the fermions for simpplicty. 

It is a standard result in open string theory 53 propagating in Abelian back- 
ground gauge Afj, and antisymmetric tensor fields that the effective target-space 
action of strings with their ends attached to a D3-brane, of interest to us here, is of 



K =A,-^A a (d f3 A^ + F^) 



(27) 




-^qO af3 F a ^r + \iq6^F a ^r Dp V 

+\mq6 a PF a pW 

-Iq^F^F^F^ + IqO^F^F^F^. 



(28) 
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Born-Infeld type: 



-BI-D3 



2^7 / d 3 x v /-Dct 4 (g^ + + a'F^) , 



(29) 



where we assume for generality that there is an induced metric on the D3 brane <? M „, 
and the D3 brane tension is determined by the Regge slope T 3 = ^ ne rea der 

should first notice that in string theory, the background fields arc the physical fields, 
connectred via the inverse Seiberg-Witten map (27) to the NC fields in this case. 
The -B-field represents the constant electric field background, and hence the source 
of NC. 

The Born-Infeld action (29) is invariant under CPT, despite the Lorentz violation 
entailed in the direction of the constant field, and this invariance was to be expected 
from the manifest validity of this symmetry in the underlying string theory. 

The expansion of the square root to first order in the NC parameter, i.e. the 
5-field in this case, reproduces the corresponding terms in (28). To see this easily, 
we may redefine a'F^ + B^ v — > F^ v . In four space-time dimensions, we have the 
identity 53 : 

-Det 4 (g^ + F^) = ^~g (l + \F^F^ - ^ {F^F^f) , ^ = \^T pa 

(30) 

Substituting F^ v = + a'F^ in (30), taking the square root, considering the 
well-known formula for the product of two n-dimensional (n = 4 in our case) anti- 
symmetric symbols 



= Det 



^ Si 1 j 1 5i 1 j 2 

^1231 &1232 



and expanding to first order in the B-field, we can easily verify the appearance of the 
two terms involving 9 a ^ oc B al3 in (28). The terms involving the dual of the Maxwell 
tensor had not been considered in ref. 51 , as they are quartic in derivatives and the 
construction of that work was to leading order in derivatives only. Neverhtless these 
terms, along with an infinity of higher derivative terms appear in the effective action 
obtained from open strings in the background of external fields. The inclusion of 
fermion backgrounds leads to the other terms in (28), to leading order in a target- 
space derivative expansion. 

By applying the reverse Seiberg-Witten map (27), then, to the fields A, ip one 
can construct in principle the NC fields and the corresponding NC action from the 
Born Infeld action, or extensions thereof 54 involving higher orders in a target-space 
derivative expansion. 
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3.3. Non-commutativity in String theory and Causality: 
photon-photon scattering and causal time delays 

In the above context of open strings in external electric field backgrounds, one 
can consider photon-photon scattering and study the effects of the induced non- 
commutativity and its connection with causality 48 . To this end, let one consider 
the Veneziano four-point amplitude, describing the scattering of two vector (spin 
one) particles, representing photon excitations, in open first-quantized string theory 
(c.f. fig. 5). In the presence of the external electric field there are non-commutative 




Fig. 5. Photon-photon scattering in open string theory (four point Veneziano amplitude). 

phases in the amplitude. The result for a world-sheet with the topology of a disc, 
where we restrict our attention here, is 47 > 48 : 

4, ~ n cfi ( K jP *(Pi A P2+P3Ap 4 ) , K l J(pi/\P4+P 3 /\P2)\ r(-2sa')r(-2ta') 

^4 ~ g s yJ^ste +^ s t e ) r(i-KW) + 

~ n eS ( K P i (PiAP2+P4Ap 3 ) , tW „i(piAp 4 +P2Ap 3 ) N \ T{-1sa')T{-1ua') , 
y s y^sut T «- SM e J r(l+2ia') " r 

n cS ( K p i (PlAp3+P2Ap 4 ) i jy-l „i(piAp 3 +p 4 Ap2)A r(-2ta') [(-W) foi\ 
~ Us y^tuC -\- J\ tu e J r(l+2sa') ' ' 

with s = 2pi -p2, t = 2pi ■ pi and u = 2p\ ■ p% the Mandelstam variables for massless 
particle states (s+t+u = 0) in the non-trivial induced space-time metric (22), which 
should be used for contraction of target-space-time indices, and we use the notation 
p A k — 8 01 (poki — koPi). The quantities K st , K' st . . . denote the standard tree- level 
kinematic factors of the amplitude, involving the photon polarization vectors . We 
thus observe that the presence of the external field induce non-commutative phases 
in the amplitude. Moreover, the perturbative amplitude expression is valid for values 
of the electric field below the critical value, as mentioned previously, since above 
that value the space-time is destabilized. 

If we look at backward scattering (u — 0) 48 , the amplitude (31) reduces to: 

A st ~ gf (K st e 2 ^ sa ' + K' st e^^ sa ') T(-2sa')T(2 S a') , (32) 

where in this case K st — ais 2 , K' st — a2S 2 , with a^, i = 1, 2 independent of s. Using 
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the identity yT(y)T(— y) = — 7r/sin(7ry), we can write: 

A ^°ff I „ „2iriEsa' , — 2-iriEsa' \ 1 ZooA 

A si - 5 S (are + a 2 e J il^W) ' ( 33 ) 

which, upon replacing s — > s + ie, e — > + , in order to treat the poles of the ampli- 
tude (at s — n/(2a') with n=integer) appropriately, and expanding the sinusoidal 
function, leads to: 

A st ~ gf [a 1 e 2m{n+ ^ sa ' + a 2 e^ l(n -^ sa '^ + 0(e) (34) 

n>0, odd 

It is important to observe that the phases are causal, in the sense that only retarded 
outgoing waves are emitted. This is to be contrasted with the non-commutative 
local field theory case, where a corresponding calculation leads to both advanced 
and retarded waves, and thus violates causality 48 . This is an important feature of 
the string calculation. 

For our purposes in this work, we are mostly interested in the fact that there is 
a series of retarded outgoing waves, with attenuating amplitudes. The first of this 
waves is emitted with a time delay of order 

A *=I=ft. ( 35 ) 

where p° denotes the total incident energy (of the two photon states in this case). 
We notice that the induced time delays survive the limit of zero electric external 
field. 

This result can be physically interpreted as follows 48 : as the two open strings in 
the diagram of fig. 5 come together, an intermediate string state is created, which 
acquires N internal oscillator excitations, growing in size from zero length to a 
maximum length L and back to zero size (oscillating N times). The whole process 
is physically possible because the intermediate string absorbs the incident energy p° 
(conservation of energy) . Since the string has a tension 1/a' one has (let us ignore 
the external field for brevity): 

P° = £ + £ (36) 

Minimizing the right-hand-side with respect to L and taking into account that the 
ends of the intermediate string state move with the speed of light in vacuo, c = 1 
(in our units), we arrive at an estimate of the time delay which is in agreement with 
the detailed (microscopic) string amplitude calculation above (in the absence of the 
external field for simplicity): 

At - a'p° . (37) 

The result (37) is consistent with the time-space uncertainty relation (20) as well 
as the standard position-momentum string corrected uncertainties 55 APAX > 
1 + a'(Ap) 2 + . . . , (in our units h = c = 1). 

An interesting question that arises at this point is whether we can measure 
experimentally such delays, and thus probe the string scale by considering photon- 
photon scattering in dense regions, such as central regions of galaxies or active 
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galactic nuclei. The reason why we need dense regions is because the causal time 
delays (35) or (37) are additive for multiple photon-photon scatterings, and hence 
the total delay is proportional to the number of scatterings a photon undergoes in 
such galactic regions before observation. Since the total delay grows proportional 
to the photon energy, the higher energy photons will be delayed more compared to 
their lower-energy counterparts, and if there is a significant number of scatterings 
in such galactic regions this may result in observable time delays of the higher 
energy photons (provided the photons are emitted "simultaneously" (within the 
experimental sensitivities)). A natural question that arises at this stage is whether 
this "stringy" cumulative delays from multiple photon-photon scatterings at the 
source of the emission of such cosmic photons could account for the observed delays 
in MAGIC 12 and/or FERMI observations 5,e . This would imply that such celestial 
sources could probe stringy uncertainties and thus play the role of "Heisenberg 
microscopes'. 

Unfortunately, for standard brane world or stringy scenarios, with conventional 
string scales a' ~ 10~ 34 m (i.e. string mass scales of order M s ~ 10 18 GeV), 
in which the intermediate space between such sources and the observation point 
is defect free, one would need unphysically high concentrations of photons in the 
source region (whose typical extent is some hundreds of Kpc (1 pc = 10 16 m)) 
in order to produce total time delays of the observed order. Even for low string 
scales of the order of TeV, the associated time delay from a string photon-photon 
scattering with total incident energy of a few TeV, is of order At ~ 10~ 27 sec, which 
again requires enormous concentrations of scattering events at the source regions in 
order to reproduce the observed high energy photon delays in MAGIC or FERMI 
observations. 

Hence, unfortunately, for standard brane world scenarios, the possibility of ce- 
lestial sources serving as stringy uncertainty probes is not realized physically. 

3.4. A stringy model for space time medium with non-trivial 
"optical" properties: D Particle "foam" 

However, as we shall argue now, such a possibility of amplification of the stringy 
uncertainty effects may be provided in cases where the bulk space between brane 
worlds is punctured by appropriate brane defects which provide the seeds for "foam" 
structures. In such a case, photons interact with (brany) space-time "defects" as 
they traverse the enormous distances, of cosmological size, between the emission 
regions and observation points on Earth or at satellites nearby. The cumulative 
(causal) time delays from such interactions, which, as we shall discuss below, assume 
the same form as the time delays of photon-photon scattering, as far as the energy 
dependence is concerned, can then become of the order of the observed delays 
in the MAGIC or FERMI Telescopes for sufficiently high concentration of space- 
time foamy defects. This is the model of the so-called "D-particle" foam, proposed 
in 17 ' 18 ' 19 ; which we now come to discuss in some detail before moving into its 
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phenomenology. In such models, the sourced regions (Active Galactic Nuclei or 
Gamma Ray Bursts) can thus play the role of " Heisenberg microscopes", probing 
stringy uncertainties but also the space-time quantum foam structure itself. As we 
shall point out, in fact, it is rather easy to falsify a large class of such models from 
the already existing experimental data in astrophysics. 




Fig. 6. A type IIA string theory model of D-particle "foam". The model consists of appropriate 
stucks of parallel D(irichlct)branes, some of which are moving in a higher-dimensional Bulk space 
time, punctured by point-like D-brane defects (D-particles). The observable Universe is repre- 
sented by one of such moving branes, compactificd to three spatial dimensions (D3-brane). As the 
D3-branc world moves, D-particles from the bulk cross it and, thus, to an observer on the brane, 
they appear as "flashing on and off" space-time foam defects ("D-particlc foam"). Photons, rep- 
resented by open strings with their ends attached on the D3 brane, interact with these defects via 
capture/recoil, and this leads to non-trivial refractive indices. The effect is therefore "classical" 
from the bulk space time point of view, but appears as an effective "quantum foam" from the 
D3-brane observer effective viewpoint. 



In 17 we have attempted to construct a brane/string- inspired model of space time 
foam which could have realistic cosmological properties. For this purpose we ex- 
ploited the modern approach to string theory 1 , involving membrane hypersurfaces 
(D(irichlet)-branes). Such structures are responsible for reconciliating (often via du- 
ality symmetries) certain string theories (like type I), which before were discarded 
as physically uninteresting, with Standard-Model phenomenology in the low-energy 
limit. 

In particular, we considered (c.f. figure 6) a ten-dimensional bulk bounded by two 
eight-dimensional orientifold planes, which contains two stacks of eight-dimensional 
branes, compactified to three spatial dimensions. Owing to special reflective prop- 
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erties, the orbifolds act as boundaries of the ninth-dimension. The bulk space is 
punctured by point- like DO-branes (D-particles), which are allowed in type IIA 
string theory (a T-dual of type I strings 56 ) we consider in 17 and here b . These are 
massive objects in string theory , with masses M s /g s , where M s is the string mass 
scale (playing the role of the quantum gravity scale in string theory), and g s < 1 
is the string coupling, assumed weak for our purposes. These objects are viewed as 
space-time defects, analogous, e.g. to cosmic strings, but these are point-like and 
electrically neutral. I have to stress at this point that, according to modern ideas 
in string theory 1 , the scale M s is in general different from the four-dimensional 
Planck-mass scale Mp = 1.2 x 10 19 GeV/c 2 , and in fact it is a free parameter in 
string theory to be constrained by experiment. The energy scale M s c 2 can be as low 
as a few TeV; it cannot be lower than this, though, since in such a case we should 
have already seen fundamental string structures experimentally, which is not the 
case. 

Supersymmctry dictates the number of D8-branes in each stack in the model, 
but does not restrict the number of DO-branes in the bulk. For dcfiniteness, we 
restrict our attention for now to the type-IIA model, in which the bulk space is 
restricted to a finite range by two appropriate stacks of D8-branes, each stack be- 
ing supplemented by an appropriate orientifold eight-plane with specific reflecting 
properties, so that the bulk space-time is effectively compactified to a finite region, 
as illustrated in Fig. 6. We then postulate that two of the D8-branes have been de- 
tached from their respective stacks, and are propagating in the bulk. As discussed 
above, the bulk region is punctured by DO-branes (D-particles), whose density may 
be inhomogeneous. When there are no relative motions of the D3-branes or D- 
particles, it was shown in 17 that the ground-state energy vanishes, as decreed by 
the supersymmetries of the configuration. Thus, such static configurations consti- 
tute appropriate ground states of string/brane theory. On the other hand, relative 
motions of the branes break target-space supersymmetry explicitly, contributing to 
the dark energy density. 

Our model assumes a collision between two branes from the original stack of 
branes (c.f. figure 6) at an early epoch of the Universe, resulting in an initial cos- 
mically catastrophic Big-Bang type event in such non-equilibrium cosmologies 57 . 
After the collision, the branes bounce back. 

It is natural to assume that, during the current (late) era of the Universe, the 
D3-brane representing our Universe is moving slowly towards the stack of branes 
from which they emanated and the configuration is evolving adiabatically. Hence 
populations of bulk D-particles cross the brane worlds and interact with the stringy 
matter on them. To an observer on the brane the space-time defects will appear 
to be "flashing" on and off. The model we are using involves eight-dimensional 
branes and so requires an appropriate compactification scheme to three spatial di- 

b One can extend the construction to phcnomcnologically realistic models of type IIB strings 19 , 
as we shall discuss in section 3.8. 
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mensions e.g. by using manifolds with non-trivial fluxes (unrelated to real magnetic 
fields) . Different coupling of fermions and bosons to such external fields breaks tar- 
get space supersymmetry, in a way independent of that induced by brane motion, 
which could be the dominant one in phenomenologically realistic models. The con- 
sequent induced mass splitting 58 ' 59 between fcrmionic and bosonic excitations on 
the brane world is proportional to the intensity of the flux field (a string generaliza- 
tion of the well-known Zeeman effect of ordinary quantum mechanics, whereby the 
presence of an external field leads to energy splittings, which are however different 
between (charged) fermions and bosons). In this way one may obtain phenomenolog- 
ically realistic mass splittings in the excitation spectrum (at TeV or higher energy 
scales) owing to supersymmetry obstruction rather than spontaneous breaking (this 
terminology which is due to E. Witten 60 , means that, although the ground state 
could still be characterised by zero vacuum energy, the masses of fermion and boson 
excitations differ and thus supersymmetry is broken at the level of the excitation- 
spectrum). The assumption of a population concentration of massive D-particle 
defects in the haloes of galaxies can lead to modified galactic dynamics 61 . Hence, 
we have an alternative scenario to standard cold dark matter, using vector insta- 
bilities arising from the splitting of strings and attachment of their free ends to a 
D-particle defect (c.f. fig. 6). 

One can calculate the vacuum energy induced on the brane world in such an 
adiabatic situation by considering its interaction with the D-particles as well as 
the other branes in the construction. This calculation was presented in 17 , where 
we refer the interested reader for further details. Here we mention only the results 
relevant for the present discussion. 

We concentrate first on D0-particlc/D8-brane interactions in the type-IIA model 
of 17 . During the late era of the Universe when the approximation of adiabatic 
motion is valid, we use a weak-string-coupling approximation g s <C 1. In such a case, 
the D-particle masses <~ M s /g s are large, i.e., these masses could be of the Planck 
size: M s /g s <~ Mp = 1.22 • 10 19 GeV or higher. In the adiabatic approximation 
for the relative motion, these interactions may be represented by a string stretched 
between the DO-particle and the D8-brane, as shown in Fig. 6. The world-sheet 
amplitude of such a string yields the appropriate potential energy between the D- 
particle and the D-brane, which in turn determines the relevant contribution to 
the vacuum energy of the brane. As is well known, parallel relative motion does 
not generate any potential, and the only non-trivial contributions to the brane 
vacuum energy come from motion transverse to the D-brane. Neglecting a velocity- 
independent term in the D0-particle/D8-brane potential that is cancelled for a D8- 
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brane in the presence of orientifold Os planes 1 c , we find: 

Ko-m = -tt£ fOT r « V^ 7 , (38) 
Vfe>8 = +8^ fOT ^»Vo7- (39) 

where u < 1 is the relative velocity between the D-particle and the brane world, 
which is assumed to be non-relativistic. We note that the sign of the effective po- 
tential changes between short distances (38) and long distances (39). We also note 
that there is a minimum distance given by: 

r min ~ vW 7 , v < 1 , (40) 

which guarantees that (38) is less than r/a', rendering the effective low-energy 
field theory well-defined. Below this minimum distance, the D0-particle/D8-brane 
string amplitude diverges when expanded in powers of (a') 2 v 2 /r . When they are 
separated from a D-brane by a distance smaller than r min , D-particles should be 
considered as lying on the D-brane world, and two D-branes separated by less than 
r m i n should be considered as coincident. 

We now consider a configuration with a moving D8-brane located at distances 
Ri{t) from the orientifold end-planes, where R\{t) + R.2{t) = Ro the fixed extent 
of the ninth bulk dimension, and the 9-density of the D-particles in the bulk is 
denoted by n*(r): see Fig. 6. The total D8-vacuum-energy density p s due to the 
relative motions is 17 : 

pZ 8 ^ D ° = - [\ n(r) *&£dr- n*(r) dr + 

/-JJi(i) 
« 4 W^* + f W n*(r)^dr + p (41) 
-is 

where the origin of the r coordinate is placed on the 8-brane world and po combines 
the contributions to the vacuum energy density from inside the band — r m ; n < r < 
r m im which include the brane tension. When the D8-brane is moving in a uniform 
bulk distribution of D-particles, we may set n* * (r) = no, a constant, and the dark 
energy density p^tai^ 011 ^ ne D8-brane is also (approximately) constant for a long 
period of time: 

Pfoti™ = -no&«(l ~v) + n^^R^tf + R 2 {tf 2a') + p Q . (42) 

Because of the adiabatic motion of the D-brane, the time dependence of Ri(t) is 
weak, so that there is only a weak time dependence of the D-brane vacuum energy 
density: it is positive if p > 0, which can be arranged by considering branes with 
positive tension. 

However, one can also consider the possibility that the D-particle density n*(r) 
is inhomogeneous, perhaps because of some prior catastrophic cosmic collision, or 

c This cancellation is crucial for obtaining an appropriate supersymmetric string ground state with 
zero ground-state energy. 
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some subsequent disturbance. If there is a region depleted by D-particles - a D-void 
- the relative importance of the terms in (42) may be changed. In such a case, 
the first term on the right-hand side of (42) may become significantly smaller than 
the term proportional to R 2 (t)/a'. As an illustration, consider for simplicity and 
concreteness a situation in which there are different densities of D-particles close 
to the D8-brane (ni oca i) and at long distances to the left and right (ni e f t , n r i g ht)- In 
this case, one obtains from (41): 

- -"local &v(l -v) + n le{t v 2 j^r (Ri(t) 2 - a') + 
nright^ 2 !^ (i? 2 (*) 2 - a') + Po ■ (43) 

for the induced energy density on the D8-brane. The first term can be significantly 
smaller in magnitude than the corresponding term in the uniform case, if the local 
density of D-particles is suppressed. Overall, the D-particle-induced energy density 
on the D-brane world increases as the brane enters a region where the D-particle 
density is depleted. This could cause the onset of an accelerating phase in the ex- 
pansion of the Universe. It is intriguing that the red-shift of GRB 090510 6 is in the 
ballpark of the redshift range where the expansion of the Universe apparently made 
a transition from deceleration to acceleration 62 . According to the above discussion, 
then, within our string foam model this may not be a coincidence 11 . 

The result (43) was derived in an oversimplified case, where the possible effects of 
other branes and orientifolds were not taken into account. However, as we now show, 
the ideas emerging from this simple example persist in more realistic structures. As 
argued in 17 , the presence of orientifolds, whose reflecting properties cause a D- 
brane on one side of the oricntifold plane to interact non-trivially with its image 
on the other side, and the appropriate stacks of D8-branes are such that the long- 
range contributions of the D-particles to the D8-brane energy density vanish. What 
remain are the short-range D0-D8 brane contributions and the contributions from 
the other D8-branes and 08 orientifold planes. The latter are proportional to the 
fourth power of the relative velocity of the moving D8-brane world: 

1) 1/ {aB, -bR 2 {t))v i (AA\ 
• / long,D8-D8„D8-08 — ^8 2 13 7r 9 a' B ' v 44 j 

where V% is the eight-brane volume, Rq is the size of the oricntifold-compactified 
ninth dimension in the arrangement shown in Fig. 6, and the numerical coefficients 
in (44) result from the relevant factors in the appropriate amplitudes of strings 
in a nine-dimensional space-time. The constants a > and b > depend on the 
number of moving D8-branes in the configuration shown in Fig. 6. If there are just 
two moving D8-branes that have collided in the past, as in the model we consider 
here, then a = 30 and b = 64 17 . The potential (44) is positive during late eras 
of the Universe as long as i? 2 (i) < 15i?o/32. One must add to (44) the negative 
contributions due to the D-particles near the D8-brane world, so the total energy 
density on the 8-brane world becomes: 
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As above, n s hort denotes the nine-dimensional bulk density of D-particles near the 
(moving) brane world. 

As in the previous oversimplified example, the transition of the D8-brane world 
from a region densely populated with D-particles to a depleted D-void causes these 
negative contributions to the total energy density to diminish, leading potentially 
to an acceleration of the expansion of the Universe. The latter lasts as long as 
the energy density remains positive and overcomes matter contributions. In the 
particular example shown in Fig. 6, R%(t) diminishes as time elapses and the D8- 
brane moves towards the right-hand stack of D-branes, so the net long-distance 
contribution to the energy density (42) (the first term) increases, tending further 
to increase the acceleration. 

It is the linear density of the D-particle defects n(z) encountered by a propa- 
gating photon that determines the amount of refraction. The density of D-particlcs 
crossing the D-brane world cannot be determined from first principles, and so may 
be regarded as a parameter in phenomenological models. The flux of D-particles is 
proportional also to the velocity v of the D8-brane in the bulk, if the relative motion 
of the population of D-particles is ignored. 

In order to make some phenomenological headway, we adopt some simplifying 
assumptions. For example, we may assume that between a redshift z < 1 and today 
(z = 0), the energy density (45) has remained approximately constant, as suggested 
by the available cosmological data. This assumption corresponds to d : 

* % = H(z)(l + z)% = - ^T^(l - v) ■ (46) 

where t denotes the Robertson- Walker time on the brane world, for which d/dt = 
—H(z)(l + z)d/dz, where z is the redshift and H(z) the Hubble parameter of the 
Universe, and we take into account the fact that dR,2{t)/dt — —v with v > 0, due to 
the motion of the brane world towards the right-hand stack of branes in the model 
illustrated in Fig. 6. Using a(t) = ao/(l + z) for the cosmic scale factor, equation 
(46) yields: 

= -H(z)(l + z)^T^(l - I ) , (47) 

where we have incorporated explicitly the speed of light in vacuo, c. Equation (47) 
may be then integrated to yield: 

n shor t(z) = n short (0) - 2T J% a * f s Jo h(z>hi+z>) whcrc e s = Vb7 . (48) 

We shall make use of this result later on, when we fit the available data from MAGIC 
and FERMI telescopes to this model. 

Before closing this section we would like to make a comment regarding the 
modification of (thermal) Dark Matter (DM) relic abundances in the foam model. As 
discussed in 63 , the quantum fluctuations of the D-particles act as sources of particle 

d We recall that, for the case of two D8 branes moving in the bulk in the model of 17 , we have 
b = 64 = 2 6 . 
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production and thus affect the respective Boltzmann equation determining the relic 
abundance of DM particles. Additional modifications to this equation are due to 
the Finsler-like 46 character of the induced space-time metric during the interaction 
of neutral particles (like the DM ones) with the foam, due to its dependence on 
D-particle recoil velocities (and hence momentum transfer) (c.f. (58) below). The 
modifications are suppressed by the square of the string scale (actually the D- 
particle mass M s /g s ) and hence for relatively high string scales (much higher than 
TeV, of interest to us here) are small and do not lead to significant constraints on 
the density of foam. 

It is important in this latter respect to make a comment regarding the nature of 
our foam: the D-particles in our approach 17 > 18 > 19 are viewed as background configu- 
rations and not as excitations of the string vacuum. In other string/brane models, 
some authors have viewed the D-particles as localized excitations of the vacuum 64 . 
In those cases, the D-particles may be considered as dark matter candidates them- 
selves and their density would be constrained by the cosmological observations on 
the DM sector to avoid overclosure of the Universe. As such they could not con- 
tribute to the refractive index. 

Indeed, for the latter property to occur, light must interact coherently with the 
D-particles, rather than scattering on them individually. Otherwise, there would 
not just be a time delay and thus an index of refraction, but the light would be 
incoherently deflected at arbitrarily large angles. Coherent scattering can only occur 
if the wavelength of light is much greater than the mean separation of the scatterers. 
In case the DO particles behave like dark matter 64 , with masses near the Planck 
scale (to account for the MAGIC delays 12 ), their number density must be less than 
2Q-20 m -3 avoid overclosure of the Universe. Thus their mean separation would 
be greater than 100 m. The gamma rays that are of interest to us here have a much 
shorter wavelength (smaller than cm) and therefore cannot experience refraction 
due to these D-particles c . 

In our D-foam model 17 ' 18 , where the DO-branes (or the compactified D3 branes 
around 3-cycles in the model of 19 ) are viewed only as background defects, such 
an issue does not arise, and the density of the foam cannot be constrained by 
overclosure of the Universe issues. As we have discussed above, in this case the D- 
particles contribute to the dark energy sector (when in motion) consistently with 
current cosmological data. The cancellation between attractive (gravitational) and 
repulsive flux forces guarantees a supersymmetric vacuum, with zero vacuum energy, 
if no relative motion of foam occurs. 



e At any rate, such super-heavy DM would have been washed out by inflation in any realistic 
cosmology, so the scenario of 64 for D-particle excitations to play the role of DM pertains to much 
lighter D-particlcs in theories with low string mass scales. 



String Quantum Gravity, LIV and Gamma-Ray Astronomy 33 



3.5. Time Delays in D-particle foam: a matter of Uncertainty 
with. . . strings attached 

After these global considerations, we can now proceed to discuss a possible origin 
of time delays induced in the arrival time of photons, emitted simultaneously from 
an astrophysical object, as a result of their propagation in the above-described D- 
particle space-time foam model. This comes by considering local interactions of 
photons, viewed as open string states in the model, with the D-particle defects. As 
we shall argue, the above model is in principle capable of reproducing photon arrival 
time delays proportional to the photon energies, of the kind observed in the MAGIC 
experiment 3 ' 12 . This microscopic phenomenon, which is essentially stringy and 
does not characterize local field theories, contributes to a sub-luminal non-trivial 
refractive index in vacuo, induced by the capture of photons or electrically neutral 
probes by the D-particle foam 18 . The capture process is described schematically in 
figure 7 and we next proceed to describe the underlying physics, which is essentially 
stringy. 
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Fig. 7. Schematic view of the capture process of an open string state, representing a photon 
propagating on a D3-brane world, by a D-particle on this world. The intermediate string state, 
indicated by thick wave lines, which is created on capture of the cnd(s) of the photon by the 
D-particle, stretches between the D-particle and the branc world, oscillates in size between and a 
maximal length of order a'p , where p° is the energy of the incident photon, and thus produces a 
series of outgoing photon waves, with attenuating amplitudes, constituting the re-emission process. 
The intermediate string state provides also the restoring force, necessary for keeping the D-particle 
roughly in its position after scattering. 



An important feature of the model is that, on account of electric charge conserva- 
tion, only electrically neutral excitations are subjected to capture by the D-particles. 
To the charged matter the D-particle foam looks transparent. This is because the 
capture process of fig. 7 entails a splitting of the open string matter state. Charged 
excitations are characterized by an electric flux flowing across the string, and when 
the latter is cut in two pieces as a result of its capture by the D-particle defect, 
the flux should go somewhere because charge is conserved. The D-particle, being 
neutral, cannot support this conservation, and hence only electrically neutral exci- 
tations, such as photons, are subjected to this splitting and the associated delays. 
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The reader should bear in mind of course that the D-particles carry other kinds of 
fluxes, unrelated to electromagnetism 1 . These are conserved separately, and it is 
for this reason that isolated D-particles cannot exist, but there must always be in 
the company of other D-branes, as in our model above, so that the relevant fluxes 
are carried by the stretched strings between the latter and the D-particles. 

For our purposes in this work it is also important to remark that the D-particles 
are treated as static when compared to photons. This is because the ends of the open 
string representing the photon move on the D3 brane world with the speed of light 
in (normal) vacuo, while the relative velocities of the D-particles with respect to 
the brane world (which propagates in the bulk space) are much lower than this. For 
instance, as discussed in 57 , to reproduce cosmological observations in this model, 
in particular the spectrum of primordial density fluctuations, which are affected by 
the relative motion of the brane world, the speed of propagation of the D3-brane 
Universe should be smaller than 10~ 4 c. 

When the end(s) of the open-string photon state are attached to the D-particle, 
there is an intermediate string state formed, stretched between the D-particle and 
the D3-brane, which absorbs the incident energy p° of the photon state to grow in 
size from zero to a maximum length L that is determined by the requirement of 
energy minimization as follows: one assumes that the intermediate string state 
needs N oscillations to achieve its maximal length L, as in the standard string- 
string scattering case, examined in section 3.3 above. Following the same logic, we 
first observe that energy conservation during the capture process of fig. 7 leads to 
relations of the form (36), from which one arrives at causal time delays for the 
outgoing photon waves of order 18 : 

At ~ a'p° . (49) 

This delay is causal and can also be obtained by considering scattering of open 
string states off D-particle backgrounds. For type IIA strings, which admit point- 
like D-particles, the situation is a straightforward extension of techniques applied 
to the open string-string scattering case of 4S , discussed in section 3.3. One as- 
sumes that the incident open string state splits into two open string states upon 
interaction with the D-particle defect, which after their scattering in the D-particle 
background, the presence of which is implemented through the appropriate Dirichlct 
world-sheet boundary conditions of the open strings, re-join to produce the outgoing 
state. The intermediate stage involves two open strings that scatter off each other. 
The (quantum) oscillations of the intermediate string state will produce a series of 
outgoing wave-packets, with attenuating amplitudes, which will correspond to the 
re-emission process of the photon after capture. The presence of the stretched string 
state, which carries the characteristic flux of the D-brane interactions, provides the 
restoring force, necessary to keep the D-particle in its position after scattering with 
the photon. 

The situation may be thought of as the stringy /brany analogue of the restoring 
force in situations in local field theories of photons propagating in media with non 
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trivial refractive indices, as discussed by Feynman 21 and reviewed above in sec- 
tion 2.1. The role of the electrons in that case (represented as harmonic oscillators) 
is played here by the D-particle defects of the space-time. The stringy situation, 
however, is more complicated, since the D-particles have an infinite number of os- 
cillatory excitations, represented by the various modes of open strings with their 
ends attached to them. Moreover, contrary to the conventional medium case, in the 
string model the refractive index is found proportional to the photon frequency, 
while the effective mass scale that suppresses the effect (49) is the quantum gravity 
(string) scale M s and not the electron mass, as in (7). The latter property can be 
understood qualitatively by the fact that the mass of these D-particle defects is of 
order 1 l/(g s \/cV) — M s /g s , where g s is the string coupling (in units of ft = c = 1). 

We remark at this stage that the above time delays are a direct consequence of 
the stringy uncertainty principles. Indeed, strings are characterized by two kinds of 
uncertainty relations: the phase-space Heisenberg uncertainty, modified by higher 
order terms in a 1 55 as a result of the existence of the minimal string length l s = \fa! 
in target space-time, 

AXAP>h+a'(AP) 2 + ... (50) 

and the time-space uncertainty relation (20). Since the momentum uncertainty 
AP < p°, we have from (50), to leading order in a' (in units h = 1): 

Ax > > w ( 51 ) 

In view of (49), we then arrive at consistency with the space-time un certainty (20), 

AX > 7& ~ S ( 52 ) 
As in the conventional string theory photon-photon scattering, reviewed in section 
3.3, these delays are causal, i.e. consistent with the fact that signals never arrive 
before they occur. Hence they are additive for multiple scatterings of photons by the 
foam defects from emission till observation. As we shall discuss in section 3.9, this 
provides the necessary amplification, so that the total delay of the more energetic 
photons can be 18 of the observed order in MAGIC and FERMI experiments. 



3.6. D-particle Recoil as External Lorentz- Violating cr-model 
Background and Non-Commutativity in D-particle foam 

In addition to this leading refractive index effect (49), there are corrections induced 
by the recoil of the D-particle itself, which contribute to space-time distortions that 
we now proceed to discuss. From a world-sheet view point, the presence of .D-particle 
recoil may be represented by adding to a fixed-point (conformal) cr-model action, 
the following deformation 42 ,43 : 

K mP = lh Ef=i fen dr «,X°9 (*°) d n X\ (53) 

where D in the sum denotes the appropriate number of spatial target-space dimen- 
sions. For a recoiling D-particle confined on a D3 brane, D = 3. 
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There is a specific type of conformal algebra, termed logarithmic conformal 
algebra 65 , that the recoil operators satisfy 42 > 43 . This algebra is the limiting case 
of world-sheet algebras that can still be classified by conformal blocks. The impulse 
operator Q(X°) is regularized so that the logarithmic conformal field theory algebra 
is respected . The conformal algebra is consistent with momentum conservation 
during recoil 42 > 43 7 which allows for the expression of the recoil velocity Ui in terms 
of momentum transfer during the scattering 

Ui=9s^, (54) 

with being the D-particle "mass" and Ap = p\ — p 2 the associated momentum 
transfer of a string state during its scattering with the D-particle. 

We next note that one can write the boundary recoil/capture operator (53) as a 
total derivative over the bulk of the world-sheet, by means of the two-dimensional 
version of Stokes theorem. Omitting from now on the explicit summation over re- 
peated i-index, which is understood to be over the spatial indices of the D3-brane 
world, we write then: 

K n ip = ihf D d 2 z ( [ Ui x°] e e &°>x*) = 



I ± 7 J D d^z(2u l )e a ^X" 



e £ {x°)+x°s £ (x°) 



d a X l (55) 



where 5 £ (X°) is an e-regularized 5-function. For relatively large times after the 
impact at X° = (which we assume for our phenomenological purposes in this 
work) , this is equivalent to a deformation describing an open string propagating in 
an antisymmetric B^ v -background corresponding to an external constant in target- 
space "electric" field, 

B M ~ Ul , By = , (X° > 0) (56) 

where the X°S(X°) terms in the argument of the electric field yield vanishing con- 
tributions in the large time limit, and hence are ignored from now on. 

To discuss the space time effects of a recoiling D-particle on an open string state 
propagating on a D3 brane world, we should consider a cr-model in the presence of 
the B-ficld (56), which leads to mixed-type boundary conditions (18) for open strings 
on the boundary dV of world-sheet surfaces with the topology of a disc. Absence 
of a recoil- velocity Ui-field leads to the usual Neumann boundary conditions, while 
the limit where — > 0, with ^ 0, leads to Dirichlet boundary conditions. 

In analogy with the standard string case in background electric fields, discussed 
in section 3.3, one obtains a non-commutative space-time if recoil of the D-particle 
is taken into account. This can be seen upon considering commutation relations 



f This can be done by using the world-sheet scale, e 2 = In (L/a) 2 , with a an Ultra- Violet scale 
and L the world-sheet area, as a regulator 42 - 43 : O e (X ) = -i Jfj^ ^j^e iu>x ° . The quantity 
s — > 0+ at the end of the calculations. 
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among the coordinates of the first quantised a- model in the background (56). As 
in the standard case of a constant electric background field, in the presence of 
a recoiling D-particle, the pertinent non commutativity is between time and the 
spatial coordinate along the direction of the recoil velocity field (for large times t 
after the impact at t = 0): 

[x 1 ,t] = i0 1 \ 0oi (= _0i O) ^ = j__j^ (57) 

where, for simplicity and concreteness, we assume recoil along the spatial X 1 di- 
rection. Thus, the induced non commutativity is consistent with the breaking of 
the Lorentz symmetry of the ground state by the D-particle recoil. The quantity 
iii = ^ and u c = ^7 is the Born-Infeld critical field. The space-time uncertainty 
relations (57) are consistent with the corresponding space-time string uncertainty 
principle (20). 

Of crucial interest in our case is the form of the induced open-string effective 
target-space-time metric. The situation parallels that of open strings in external 
electric field backgrounds, discussed in refs. 48 ' 47 J and reviewed in section 3.1. Hence, 
the effective open-string metric, gr°P en . electric 5 which is now due to the presence of 
the recoil-velocity field u, whose direction breaks target-space Lorentz invariance, 
is obtained by extending appropriately the result (22) to the background (56): 

ff opc„,olcctric = (1 _ -2) ^ ; ^ = 0,1 

ff opc„,olcctric = r] ^ li y= al j Qther valueg _ ( 5g ) 

For concreteness and simplicity, we considered a frame of reference where the matter 
particle has momentum only across the spatial direction X 1 , i.e. 7^ p\ = p | 
u i 1V1 = P3 = 0. Moreover, as in the standard case of strings in an electric field 
background, there is a modified effective string coupling 47,48 (c.f. (23)): 

gf = : g s (l- u 2 ) 1/2 (59) 

The fact that the metric in our recoil case depends on momentum transfer variables 
implies that D-particle recoil induces Finsler-type metrics 46 , i.e. metric functions 
that depend on phase-space coordinates of the matter (photon) state. 

As already mentioned, the presence of the critical "background field" u c is as- 
sociated with the destabilization of the vacuum when the field intensity approaches 
the critical value. Since in our D-particle foam case, the role of the 'electric' field is 
played by the recoil velocity of the D-particle defect, the critical field corresponds 
to the relativistic speed of light, in accordance with special relativistic kinematics, 
which is respected in string theory, by construction. On account of (54), then, this 
implies an upper bound on the induced momentum transfer, and hence on the avail- 
able momenta, for the effective field theory limit to be valid. Indeed, if we represent 
Ap in (54) as a fraction of the incident momentum p\, Ap = rpi, r < 1, then 
the condition that the recoil velocity of the D-particle is below the speed of light 
in vacuo, as required by the underlying consistency of strings with the relativity 
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principle, implies 22 

Ap=r Pl < %=>Pi < % • (60) 

When the incident momentum approaches the order of this cutoff, the effective string 
coupling (59) vanishes, while above that value the coupling becomes imaginary, 
indicating complete absorption of the string state by the D-particle. The space- 
time distortion due to recoil is so strong in such a case that there is no possibility 
of re-emergence of the string state, the defect behaves like a black hole, capturing 
permanently the string state. 

It is important to notice that in modern string theory the quantities M s ,g 8 are 
completely phenomcnological. In fact, it is possible to construct phenomenologically 
realistic string/brane-Universe models, in the sense of being capable to incorporate 
the Standard Model particles at low energies, with string scales in such a way 
that M s /g s is significantly lower than the Planck scale. For instance, there are 
constructions 92 in the large extra dimension framework, for which M s /g s is of the 
order of 10 19 eV, i.e. of order of the conventional GZK cutoff in Lorentz invariant 
particle physics models. Thus, by embedding the D-particle foam to such models, 
one may have the appearance of a novel type Gretisen Zatsepin Kusmin (GZK) 
cutoff 66 , of similar order to the one obtained from conventional Lorentz invariance 
arguments, but of quite different origin: here it is the subluminal nature of the recoil 
velocity of the foam that sets the new upper bound in momentum transfer. This is 
not fine-tuning in our opinion, but indicates the appearance of a new upper bound 
in momenta, related implicity to the underlying Lorentz invariance of the string 
theory, which is broken spontaneously by the recoiling D-particle background. This 
point will become of importance later on, when we discuss constraints of our model 
coming from ultra- high-energy/infrared photon/photon scattering 20 . 

Before closing this subsection we make an important remark. The induced metric 
(58) will affect the dispersion relations of the photon state: 

P »p v g°l en ' electlic = . (61) 

However, because the corrections on the recoil velocity Ui are quadratic, such modi- 
fications will be suppressed by the square of the D-particle mass scale. On the other 
hand, the presence of the D-particle recoil velocity will affect the induced time de- 
lays (49) by higher-order corrections of the form, as follows by direct analogy of our 
case with that of open strings in a constant electric field 48 : 

^^with D— foam recoil velocity — ex • (62) 

Thus, the D-particle recoil effects are quadratically suppressed by the D-particle 
mass scales, since (c.f. (54)) vn oc g s Api/M s , with Api the relevant string-state 
momentum transfer. However, the leading order delay effect (62), obtained formally 
by considering the limit of vanishing recoil velocity, is linearly suppressed by the 
string scale, and thus the induced time delays are disentangled in the string foam 
model from the modified Finslcr dispersion relations (61). This is important for the 
phenomenology of the model, as we shall discuss below. 
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3.7. Lack of a Local Effective Field Theory Formalism in recoiling 
D-particle Foam 

An important comment arises at this stage concerning the construction of a possible 
local effective field theory action in the case of D-particle foam, when recoil of the 
D-particle defect is considered. In view of the formal analogy of the problem with 
that of an external electric background field, one is tempted to apply the same 
considerations as those in section 3.2 leading to the non-commutative effective action 
(28) for the case of quantum electrodynamics in the presence of D-foam background, 
with the 8 a P parameter being replaced by 6° l = U{ = g s ^r, the D-particle recoil 
velocity. 

However, this is not correct. As a result of the momentum-transfer dependent 
nature of the non-commutativity parameter in the D-foam case, it is not possible to 
write down the effective action in target space as a power series of local quantum 
operators. The momentum transfer is not represented by such operators when taking 
the Fourier transform. 

Thus writing down a local effective action for the recoiling D-particle in inter- 
action with, say, a photon, is not possible. This has important consequences for 
phenomenology. Indeed, the time delays (62), which are associated with the stringy 
uncertainty, are found proportional to the incident energy of the (split) photon state 
and are thus linearly suppressed by the string scale. The associated refractive index 
is therefore linearly suppressed, but this cannot be interpreted as an average prop- 
agation of photons in the context of some local effective action. As we have seen 
above, the associated anomalous photon dispersion terms in this case, induced by 
the Finsler metric (58), are quadratically suppressed by the string scale. 

Therefore, any analysis, such as those involving high energy cosmic rays 20 using 
an effective linear dispersion relation, obtained from a local effective action, does 
not apply to our problem. Moreover, it is a generic feature of any local effective 
theory, that is a theory on flat space times involving higher derivatives operators, to 
yield birefringence, since the corresponding modified dispersion relations for photons 
reduce at the end of the day to a solution of the photon frequency as a function of 
the wave number k which is obtained from a quadratic equation. The latter admits 
two solutions with different propagation for the two photon polarizations. 

This is not the case of the uncertainty related time delay (62), which as we have 
seen is independent of the photon polarization 18 . 

3.8. D-particle "foam" in Type LIB strings 

The discussions in the previous subsections referred formally to D-brane defects of 
zero spatial size (point-like D-particles), which are allowed in type IIA string theory. 
One can do a similar analysis for (the phenomenologically more realistic, from the 
point of view of incorporating Standard Model Physics in the low-energy limit) 
type IIB strings 19 , which, though, do not contain D-particles (DO branes) . In such 
a case the "foamy" defects can be constructed by wrapping up D3 branes around 
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three cycles. The role of brane worlds in such type IIB string models is played by 
appropriately compactified D7 branes. There are some technical subtleties in this 
case which are worthy of discussing, and which lead to some important differences 
from the type IIA D-particle foam model. 

Let us consider the Type IIB string theory with D3-branes and D7-branes where 
the D3-branes are inside the D7-branes. The D3-branes wrap a three-cycle, and 
the D7-branes wrap a four-cycle. Thus, the D3-branes can be considered as point 
particles in the Universe, i.e., the D-particles, while the Standard Model (SM) 
particles are on the world- volume of the D7-branes. 




(c) 



Fig. 8. (a): the splitting/capture/re-emission process of a (generic) matter string by a D-particle 
from a target-space point of view, in the type IIB D-foam model, (b) The same process but for 
photons (or in general particles in the Cartan subalgebra of the Gauge group in this (intersecting) 
brane world scenario, (c): The four-point string scattering amplitude (corresponding to the parts 
inside the dashed box of (a)) between the constituent open ND strings of the splitting process. 
Latin indices at the end-points of the open string refer to the brane worlds these strings are 
attached to. 



For simplicity, we assume that the three internal space dimensions, which the 
D3-brane wraps around, are cycles S 1 x S 1 x S 1 , and we denote the radius of the 
i— th cycle S 1 as Ri. The mass of the D3-brane is 67 

,, R1R2R3 , aQ . 

D3 = af A > ( 63 ) 

whereby g s is the string coupling, and l s is the string length, i.e., the square root 
of the Regge slope £ s = \[a! . If Ri < £ s , we can perform a T-duality transformation 
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along the i— th S 1 as follows 



R ^f.> (64) 



and then we obtain i?^ > £ s . In string theories with compact internal space di- 
mension^), there are Kaluza-Klein (KK) modes and string winding modes. Under 
T-duality, the KK modes and winding modes are interchanged. In particular, these 
two theories are physically identical 67 . Therefore, without loss of generality, we can 
assume Ri > l s . Choosing Ri = 10£ s and g s <~ 0.5, we obtain Md3 ~ 2000/4, and 
then the mass of the D3-brane can be much larger than the string scale. Therefore, 
anticipating the result (62), we shall ignore to a first approximation the recoil of 
the D-particles. Their inclusion (as small, perturbative corrections) is straightfor- 
ward, according to the discussion in the previous section, and we shall come back 
briefly to this issue at the end of the article. In this setting, the particles (called 
ND particles) arising from the open strings between the D7-branes and D3-branes 
which satisfy the Neumann (N) and Dirichlet (D) boundary conditions respectively 
on the D7-branes and D3-branes, have gauge couplings (with the gauge fields on 
the D7-branes) of the form: 

4- = %, (65) 

where g-j are the gauge couplings on the D7-branes, and V denotes the volume of 
the extra four spatial dimensions of the D7 branes transverse to the D3-branes 68 . 
Because the Minkowski space dimensions are non-compact, V is infinity and then 
537 is zero. Thus, the SM particles have no interactions with the ND particles on 
the D3-brane or D-particle. 

To have non-trivial interactions between the particles on D7-branes and the 
ND particles, we consider a D3-brane foam, i.e., the D3-branes are distributed 
uniformly in the Universe, as a background. We assume that the U43 is the average 
three-dimensional volume that has a D3-brane in the Minkowski space dimensions, 
and R' is the radius for the fourth space dimension transverse to the D3-branes. 
In addition, in the conformal field theory description, a D-brane is an object with 
a well defined position. In contrast, in the string field theory, a D-brane is a "fat" 
object with thickness of the order of the string scale. In particular, the widths of 
the D-brane along the transverse dimensions are about 1.55^ s , as follows from an 
analysis of the tachyonic lump solution in the string field theory which may be 
considered as a D-brane 69 . Thus, our ansatz for the gauge couplings between the 
gauge fields on the D7-branes and the ND particles is 19 : 

1 _ Va 3 R' it _ VasR' 1 (ac\ 

We denote a generic SM particle as an open string ab with both ends on the 
D7-branes (c.f. fig. 8). For a = b, we obtain the gauge fields related to the Cartan 
subalgebras of the SM gauge groups, and their supersymmetric partners (gauginos), 
for example, the photon, Z° gauge boson and the gluons associated with the A 3 and 
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A 8 Gell-Mann matrices of the SU(3)c group. For a ^ b, we obtain the other parti- 
cles, for example, the electron, neutrinos, Wrf boson, etc. As indicated in fig. 8(a), 
when the open string ab passes through the D3-brane, it splits and becomes two 
open strings (corresponding to the ND particles) ac and cb with one end on the 
D7-brane (a or b) and one end on the D3-brane (c or c). Then, we can have the 
two-to-two process and have two outgoing particles arising from the open strings 
ac and cb. Finally, we can have an outgoing particle denoted as open string ab. In 
particular, for a — b, at leading order we can have the s-channel process depicted 
in fig. 8(b). 

The time delays arise from the two-to-two process in the box of fig. 8(a), with 
the corresponding string diagram give in fig. 8(c). To calculate the time delays, 
we consider the four-fermion scattering amplitude and use the results of rcf. 70 for 
simplicity. We can discuss the other scattering amplitudes similarly, for example, 
the four-scalar scattering amplitude, the results for the delay are the same. The 
total four-fermion scattering amplitude is obtained by summing up the various 
orderings 70 : 

Aotai = A(l, 2, 3, 4) + .4(1, 3, 2, 4) + .4(1, 2, 4, 3) , (67) 

where 

.4(1, 2, 3, 4) = A(l, 2, 3, 4) + A(4, 3, 2, 1) . (68) 
A(l,2,3,4) is the standard four-point ordered scattering amplitude 
(27r) 4 <5( 4 )(^fc a )A(l,2,3,4) = ^/> 

a 

(v (1) (0,fc 1 )V (2) (a;,fc 2 )V (3) (l,fc3)V (4) (oo,A : 4)) , (69) 

where ki are the space-time momenta, and we used the SL(2,R) symmetry to fix 
three out of the four Xi positions on the boundary of the upper half plane, represent- 
ing the insertions of the open string vertex fermionic ND operators , i = 1, . . . 4, 
describing the emission of a massless fermion originating from a string stretched 
between the D7 brane and the D3 brane 70 . 

The amplitudes depend on kinematical invariants expressible in terms of the 
Mandclstam variables: s = — (k\ + k 2 ) 2 , t = — (fc 2 + k 3 ) 2 and u = — (fci + k 3 ) 2 , 
for which s + t + u = for massless particles. The ordered four-point amplitude 
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.4(1,2,3,4) is given by 

-g s l 2 s fdxx- 1 -'* (1-.T)- 1 -"' 
Jo 

V 1 »7 ( ,u (2) « (4 ytt (3) (l - + « (1) 7^ (4) « (2) 7 M « (3) ^ 

jA E e_7rr " 2 R ' 2/ ^> ■ ( 7 °) 

nez 

where F(x) = F(l/2; 1/2; 1; x) is the hypergeometric function, t(x) — F(l — 
x)/F(x), ji and Jj with i = 1, 2, 3, 4 are indices on the D7-branes and D3-branes, 
respectively, and is 

(1-554) 4 rrn 
= ~V^B r ' (71) 

in the notation of 70 , u is a fermion polarization spinor, and the dependence of the 
appropriate Chan-Paton factors has been made explicit. Thus, taking F(x) ~ 1 wc 
obtain 

4(1,2,3,4) oc g s t\ {ttlu^h^u^h^ 

4(1,3,2,4) oc ^ (^(D 7aiU (3)^(4) 7 m u (2) 
+uilu^u^ 3 h^ ) ) x ^gfe^ , 
4(1,2,4,3) oc g s l\ (u^(i) 7AiW (2)y(3) 7 M u (4) 

+^Sm« (3) * V« (4) ) x ^ggg^ , (72) 

where the proportionality symbols incorporate Kaluza-Klcin or winding mode con- 
tributions, which do not contribute to the time delays. Technically, it should be 
noted that the novelty of our results above, as compared with those of 70 , lies on 
the specific compactification procedure we adopted, and the existence of a uniformly 
distributed population of D-particlcs (foam), leading to (66). It is this feature that 
leads to the replacement of the simple string coupling g s in the Veneziano amplitude 
by the effective D3-D7 effective coupling gf 7 s . 

Similarly to the discussion in Ref. 48 , time delays arise from the amplitude 
4(1, 2, 3, 4) by considering backward scattering u = 0. Noting that s + 1 + u = for 



g We stress once again that this is only part of the process. The initial splitting of a photon into 
two open string states and their subsequent re-joining to form the re-emitted photon both depend 
on the couplings <?| 7 , so there are extra factors proportional to n in a complete treatment, which 
we do not discuss here, as they are not relevant for the estimation of the time delays. 
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massless particles, the first term in .4(1, 2, 3, 4) in Eq. (72) for u — is proportional 
to 

te s T{- s e s )T{-til) = -se 2 s T(-sel)T( S £l) 

= sm(irs£ 2 s ) ' ( 73 ) 

It has poles at s — n/l 2 , with n an integer. The divergence of the amplitude at the 
poles is an essential physical feature of the amplitude, a resonance corresponding to 
the propagation of an intermediate string state over long space-time distances. To 
define the poles we use the correct e prescription replacing s — > s + ie, which shift 
the poles off the real axis. Thus, the functions l/sin(7rs^) can be expanded as a 
power series in y: 

y = e ™<~ e . (74) 

On noting that s = E 2 , where E is the total energy of the incident string state (c.f. 
fig. 8), we obtain the time delay at the leading order 

At = E£ 2 S . (75) 

Let us discuss the time delays for concrete particles. We will assume that r\ is 
a small number (c.f. below). Then, for the gauge fields (and their corresponding 
gauginos) which are related to the Cartan subalgebras of the SM gauge groups, 
all the amplitudes 4(1, 2, 3, 4), 4(1, 3, 2, 4), and 4(1, 2, 4, 3) will give the dominant 
contributions to the total amplitude due to j\ = ji . Thus, they will have time 
delays as given in Eq. (75). The resulting delay for photon is independent of its 
polarization, and thus there is no birefringence, thereby leading to the evasion of 
the relevant stringent astrophysical constraints 71 > 72 > 73 j as we shall see later on. 

However, for the other particles, we have j\ ^ j2, and then only the amplitude 
4(1, 3, 2, 4) gives dominant contribution. Considering backward scattering 48 u = 
and s + t + u = 0, we obtain 

4(1,3,2,4) oc g s i 2 s (^^iWVu^ 

__L^2-(i) 7 ^(4)-( 3 ) 7 ^(2)^ _ (?6) 

Because they are just the pole terms, we do not have time delays for other particles 
with ji ^ j2 at the leading order, for example, boson, electron, and neutrinos, 
etc. At order r] (0(7/)), we have time delays for these particles, which arise from the 
forth line in Eq. (70). 

These time delays are suppressed by terms of order r\ compared to the corre- 
sponding time delays for photons or particles in the Cartan subalgebra of the gauge 
group of the model. For some characteristic values of this parameter we note that, 
quite easily, and within natural ranges of the parameters of the model, T] can be of 
order 

V < !0~ 6 • (77) 
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For example, by taking the following Va3 and R' in Eq. (71) 

Vas ~ (104) 3 , R' ~ 3384 • (78) 

one achieves the above range. 

Although, as we have already mentioned, the time delays calculated from string 
amplitudes in the above way do not necessarily correspond to modified dispersion 
relations for the probe, nevertheless we note that if one attempts to interpret the 
associated time delays for the electron in this model as modified dispersion relations, 
one would write E 2 = p 2 + m\ — rip 3 /M st , with Mg t an effective quantum gravity 
scale. For Mst of the order of a string scale in the range of 10 18 GeV, as required 
by the MAGIC fit of the observed photon delays 12 , and ij in the range (77), the 
constraints on the electron dispersion from the Crab Nebula synchrotron radiation 
observations 74 ' 75 would be easily satisfied. However, as already mentioned, the 
ordered amplitude (70) describes only part of the process, and there are extra factors 
of rj coming from the initial splitting of the photon into two open string states. In 
such a case, the scale Mst would be inversely proportional to them, and so the 
constraints from Crab Nebula can be satisfied for much larger values of r\ than 
(77). Hence, as in the case of type IIA string foam, synchrotron radiation does not 
provide stringent constraints for type IIB D-foam. 

3.9. Multiple Photon- D -particle Scatterings in the Foam and Total 
(observed) Time Delays 

The above-discussed time delays (49) (or (75)) pertain to a single encounter of a 
photon with a D-particle. In case of a foam, with a linear density of defects n* /y/o7 ', 
i.e. n* defects per string length, the overall delay encountered in the propagation 
of the photon from the source to observation, corresponding to a traversed distance 
D, is: 

At to ^ = a'p°n*^ = ^n*D , (79) 

where p° denotes an average photon energy. When the Universe's expansion is taken 
into account, one has to consider the appropriate red-shift- z dependent stretching 
factors which affect the measured delay in the propagation of two photons with 
different energies as well as the Hubble expansion rate H{z) 13 ' 14 . Specifically, in a 
Roberston- Walker cosmology the delay due to any single scattering event is affected 
by: (i) a time dilation factor 14 (1+z) and (ii) the redshifting 13 of the photon energy 
which implies that the observed energy of a photon with initial energy E is reduced 
to E ohs = E /(l + z). Thus, the observed delay associated with (49) is 13 > 14 : 

St ohs = (1 + z)5t Q = (1 + z) 2 Vc7E ohs . (80) 

For a line density of D-particles n(z) at redshift z, we have n(z)d£ — n(z)dt defects 
per co-moving length, where dt denotes the infinitesimal Robertson- Walker time 
interval of a co-moving observer. Hence, the total delay of an energetic photon in a 
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co-moving time interval dt is given by n(z)(l + z) 2 C\/a'E b s dt. The time interval 
dt is related to the Hubble rate H(z) in the standard way: dt = — dz/[(l + z)H(z)]. 
Thus, from (80) we obtain a total delay in the arrival times of photons with energy 
difference AE, which has the form considered in 13 ' 15 7 namely it is proportional to 
AE and is suppressed linearly by the quantum gravity (string) scale, M s : 

(At) ohs = ^H- 1 Cn*{z) - {1+z)dz (81) 

V Jobs Ms o JO V ; ^/0 A +O m (l+ z )3 V > 

where z is the red-shift, H is the (current-era) Hubble expansion rate, and we 
have assumed for concreteness the ACDM standard model of cosmology, with fij = 
representing the present-epoch energy densities (pi) of matter (including dark 
matter), Ct m , and dark vacuum energy, Q,\, in units of the critical density p c = 
g^Q N of the Universe (Gat is the Newton's gravitational constant), that is the 
density required so that the Universe is spatially fiat. The current astrophysical 
measurements of the acceleration of the Universe are all consistent with a non 
zero Cosmological-Constant Universe with Cold-Dark-Mattcr (ACDM Model), with 
A ~ 73% and Q m ~ 27%. 

Notice in (81) that the essentially stringy nature of the delay implies that the 
characteristic suppression scale is the string scale M s , which plays the role of the 
quantum gravity scale in this case. The scale M s is a free parameter in the modern 
version of string theory, and thus it can be constrained by experiment. As we have 
discussed in this article, the observations of delays of energetic (TeV) photons from 
AGN by the MAGIC telescope 3 can provide such an experimental way of constrain- 
ing n* /M s in (81). For AE <~ 10 TcV, for instance, the delay (81) can lead to the 
observed one of order of minutes, provided M s /n* ~ 10 18 GeV (in natural units 
with c = 1) 12 . This implies natural values for both n* and M s , although it must 
be noted that n* is another free parameter of the bulk string cosmology model of 
17 , considered here. In general, n*(z) is affected by the expansion of the Universe, 
as it is diluted by it. This depends on the bulk model and the interactions among 
the D-particles themselves. For redshifts of relevance to the MAGIC experiments, 
z = 0.034 <C 1, one may ignore the ^-dependence of n* to a good approximation. 

The total delay (81) may be thought of as implying 21 an effective subluminal 
refractive index n(E) of light propagating in this space time, since one may assume 
that the delay is equivalent to light being slowed down due to the medium effects. 
On account of the theoretical uncertainties in the source mechanism, however, the 
result of the AGN Mkn 501 observations of the MAGIC Telescope translate to 
upper bounds for the quantity n* /M s in (81), which determines the strength of the 
anomalous photon dispersion in the string/D-particle foam model. 

In view of the above discussion, if the time delays observed by MAGIC can 
finally be attributed partly or wholly to this type of stringy space-time foam, and 
therefore to the stringy uncertainties, then the AGN Mkn 501, and other such 
celestial sources of very high energy photons, may be viewed as playing the role of 
Heisenberg microscopes and amplifiers for the stringy space-time foam effects. 
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We next proceed to discuss the phenomenology of the D-foam model, and see 
whether it can survive the stringent constraints on Lorentz Violation imposed by 
a plethora of astrophysical measurements. To understand better the situation, and 
demonstrate why the exotic interpretation put forward in this article could be a 
viable hypothesis, we review briefly first the conventional Astrophysics mechanisms 
for cosmic particle acceleration, for which currently there is no consensus among 
the astrophysicists. 



4. High-Energy Gamma Ray Astronomy &; Fundamental Physics 

In this part of the review we shall discuss experimental tests of Lorentz symmetry 
using high energy cosmic probes. Then we shall try to see what these observations 
imply for the Lorentz- violating model of stringy (D-particle) space-time foam dis- 
cussed in previous sections and how the latter compare with other Lorentz violating 
field theories available to date. The important data for our purposes are delayed 
arrivals of high energy photons (compared to lower energy ones) from distant cosmic 
sources. In particular, currently there are the following "anomalous" photon events: 

• MAGIC Telescope (source: Active Galactic Nucleus (AGN) Mkn 501, red- 
shift z=0.034), Highest Energy 1.2 TeV Photons, Observed Delays of 
O(min) 3 < 12 

• HESS Telescope (source AGN PKS 2155-304, redshift z=0.116), Highest 
Energy 10 TeV photons, Originally claim no observed time lags 10 . 

• FERMI Satellite (source Gamma Ray Burst (GRB) 090816C, redshift 
z=4.35), Highest Energy Photon 13.2 GeV, 4.5 s time-lag between E > 
100 MeV and E < 100 KeV Observed Time Delay 16.5 sec 5 . 

• FERMI Satellite (source GRB 090510, redshift z=0.9), Highest Energy Pho- 
ton 31 GeV, several photons in range 1-10 GeV 6 ; this is a short and intense 
GRB. Observed Time Delays < 1 sec h . We also notice, out of curiosity, that 
the location of this GRB is in the ball park of the ranges of redshift where 
the cosmic deceleration/acceleration transition for our Universe occurs. 

• FERMI Satellite (source GRB 09092B, redshift z=1.822), Highest Energy 
Photon 33.4 GeV, Observed Time Delay 82 sec after GMB trigger, 50 sec 
after end of emission 7 . 

To understand in detail the possible explanations of the above observations, let 
us first recapitulate the up to date knowledge on the production of very high energy 
gamma rays in the Universe 76 and ways of cosmic acceleration. I must stress once 
again that there is no consensus among the astrophysical community as regards 
the various ways of particle acceleration at various regions in the Universe, and 



h However, in our opinion there are some ambiguities on the exact time emission of the highest 
energy photons, since the measurements on the precursor to the GRB are not yet conclusive as far 
as we understand. The assumption in this review, though, is that this measurement stands and 
we try to interpret it using our model in conjunction with the rest of the data. 
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in fact it is most likely that there are several mechanisms taking place, depend- 
ing on the source. The MAGIC and FERMI observations added to this puzzle, in 
particular why the more energetic photons from the respective cosmic sources ar- 
rived later than the lower energy ones with the observed time delays. Is there any 
correlation between the time delays and the photon energies ? Our understanding 
of the source mechanism for the production of high energy particles from cosmic 
sources will improve only by making more and more precision measurements from 
a variety of celestial sources, a process currently under way through the operation 
(or construction) of several terrestrial and extra-terrestrial facilities of high-energy 
astrophysics. 

It is these uncertainties in the conventional astrophysics of the sources that allow 
for speculations that fundamental physics, such as photon propagation in a quantum 
gravity "medium", might play a significant role on the MAGIC effect. However, 
if the latter has a chance of being true it must respect all the other stringent 
astrophysical constraints on Lorentz invariance that are currently available. As we 
have already discussed in section 3.4, it seems that, at present, only a specific 
string theory model of quantum foam, in which only photons are not transparent 
to the foam effects, could stand up to this chance. We shall come back to a detailed 
phenomenology of the model, later on. At present, we proceed to a brief review of 
the basic astrophysical mechanisms for cosmic particle acceleration available today 
in order to appreciate the difficulties in using them as possible explanations of the 
above-mentioned observed time delays of the more energetic photons. 

4.1. Conventional Astrophysics mechanisms for cosmic 
Very-High-Energy (VHE) Gamma-Ray production 

Gamma rays constitute the most interesting part of the spectrum of active galactic 
nuclei (AGN). An AGNs is a compact region at the centre of a galaxy, with much 
higher than normal luminosity over some or all of the electromagnetic spectrum (in 
the radio, infrared, optical, ultra-violet, X-ray and/or gamma ray wavebands). The 
radiation from AGN is believed to be a result of accretion on to the supermassive 
black hole at the centre of the host galaxy. AGN are the most luminous persistent 
sources of electromagnetic radiation in the universe, and as such can be used as a 
means of discovering distant objects; their evolution as a function of cosmic time 
also provides constraints on cosmological models. 

Gamma Rays with energies higher than 20 MeV and up to TeV have been 
observed today from such AGNs. It is customary (although somewhat arbitrary) to 
classify these Gamma Rays according to their energies as follows: 

• (i) High-Energy Gamma rays: with energies from 20 MeV - 100 GeV , 

• (ii) Very High-Energy Gamma rays: with energies from 100 GeV - 30 TeV , 

• (iii) Ultra High-Energy Gamma rays: with energies from 30 GeV - 30 PeV , 

• (iv) Extremely High-Energy Gamma rays: with energies from 30 PeV - ?. 
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Theoretically, the last category incorporate energies up to the ultraviolet cutoff 
energy scale (Planck-scale energies 10 19 GeV) that defines the structure of low- 
energy field theories as we know them. 

The production of very high energy gamma rays is still not understood well, and 
constitutes a forefront of research on galactic and/or extragalactic astrophysics. 




Fig. 9. The basics of a cosmic accelerator model. Very high energy Photon production from Grav- 
itational Energy conversion (Penrose process) in AGNs, believed to take place in AGN Mkn501. 
Relativistic matter, such as electrons, are beam ejected as a result of the enormous gravitational 
energy available during the collapse process forming a black hole at the center of the galaxy. Such 
electrons, undergo subsequent synchrotron radiation due to their interaction with the magnetic 
fields existing in the galactic region, and eventually inverse Compton scattering (IC), i.e. inter- 
actions of these very high energy electrons with low-energy photons (say of eV-keV energies), to 
produce TeV photons from Mkn501 observed by MAGIC (see fig. 10). This combined process is 
called Synchrotron-self-Compton process (Figure taken from Ref. 77 ). 



At present there are three major categories believed responsible for the produc- 
tion of very- and ultra- high-energy Gamma rays: 



• (i) Photons from conversion of gravitational energy in Active Galactic Nu- 
clei (AGN) (c.f. fig. 9), 

• (ii) From self-annihilation of Dark Matter and 

• (iii) From decays of exotic massive particles (with masses of order 10 15 — 10 16 
GeV/c 2 ), appearing in Grand Unifying Models Beyond the Standard Model 
(such as string-theory inspired models), in the very Early Universe. 
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Fig. 10. One of the suggested models (leptonic acceleration) for the production of very high energy 
gamma rays: Synchrotron-self-Compton (SSC), involving synchrotron (Sy) and inverse Compton 
(IC) scattering for the same electron, provides a mechanism for the production of photons with 
energies up to several TeV. In the middle figure the typical energy spectrum is indicated, with 
the characteristic double peak of the SSC mechanism, one of the synchrotron radiation and the 
other for the IC scattering. The Inverse Compton spectrum has a peak at TeV energies, as those 
observed by MAGIC when looking at the AGN Mkn 501 (lower picture, from Ref. 78 ). 
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In what follows we shall restrict ourselves to (i), which most likely pertains to 
the production of very high energy Gamma rays observed in AGN such as Mkn 501 
observed by the MAGIC Telescope. It is believed today that the centres of galaxies 
contain massive black holes due to matter collapse 79 , with typical masses in the 
range 10 6 — 10 9 solar masses. AGN therefore are celestial systems with very high 
mass density, and it has long been assumed that they consist of a massive black 
hole, of say 10 8 solar masses or more, accreting the gas and dust at the center of 
the galaxy. The gravitational energy liberated during accretion onto a black hole is 
10% of the rest mass energy of that matter and is the most efficient mass-energy 
conversion process known involving normal matter. Collapsing matter towards this 
massive central galactic object releases gravitational energy (Gravitational Energy 
Conversion) and results in spectacular relativistic material jet emissions. Since the 
accreting matter has in general a non-trivial angular momentum, angular momen- 
tum conservation is responsible for the matter orbiting the black hole and, through 
energy dissipation, the formation of a material (flat) accretion disk. This also re- 
sults in material jets with ultra-relativistic particles outflowing the accretion plane 
(fig. 9). 

In addition, since the black holes at the centres of the AGN are probably rotating 
(as a result of having a non-zero angular momentum (Kerr type) 79 , due to angular 
momentum conservation in the formation (collapse) process), one might also spec- 
ulate that part of the relativistic jet might be due to the so-called Penrose process, 
which allows the extraction of energy from a rotating black hole 80 . The extraction 
of energy from a rotating black hole is made possible by the existence of a region of 
the Kerr spacetime called the ergoregion, in which a particle is necessarily propelled 
in locomotive concurrence with the rotating spacetime. In the process, a lump of 
matter enters into the ergoregion of the black hole and splits into two pieces, the 
momentum of which can be arranged so that one piece escapes to infinity, whilst 
the other falls past the outer event horizon into the hole (a rotating black hole has 
two event horizons, an outer and an inner one). The escaping piece of matter can 
possibly have greater mass-energy than the original infalling piece of matter. In 
summary, the process results in a decrease in the angular momentum of the black 
hole, leading to a transference of energy, whereby the momentum lost is converted 
to energy extracted. The process obeys the laws of black hole mechanics. A conse- 
quence of these laws is that if the process is performed repeatedly, the black hole 
can eventually lose all of its angular momentum, becoming rotationally stationary. 

The particles in the relativistic jet undergo acceleration, but currently the per- 
tinent mechanism is a matter of debate and active research. Most likely it depends 
on the source. In general there are two generic ways of cosmic acceleration. 

Conventional Astrophysics Mechanisms for Cosmic Acceleration 

Leptonic Acceleration: Among the particles in the jet are charged electrons, 
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whose paths are curved as a result of the existing magnetic fields in the galactic 
regions, which accelerate the electrons. The curved path of a charged object im- 
plies synchrotron radiation, as a result of energy conservation. Moreover, in AGN's 
like Mkn 501, the same (high-energy) electron can also undergo inverse-Compton 
scattering with low-energy photons (with energies of order eV, e.g. photons of the 
cosmic microwave background radiation that populate the Universe as remnants of 
the Big-Bang). The terminology inverse-Compton (IC) scattering refers here to the 
fact that, contrary to the conventional Compton photon-electron scattering, here 
it is the electron which is the high energy particle, and whose loss of energy is 
converted to outgoing radiation (c.f. fig. 10). This IC outgoing radiation can have 
very high energies. In fact, the Compton spectrum peak can be at several TeV ener- 
gies, as observed for the AGN Mkn 501. This combined process, whereby the same 
electron that is responsible for synchrotron radiation in AGN also undergoes IC 
scattering to produce high energy photons is known as Synchrotron-self-Compton 
(SSC) mechanism 81 , and is believed 76 -with some variations to be discussed below 
- that is responsible for the production of the very high energy photons observed in 
the AGN Mkn 501 (c.f. fig. 10). This is the so-called Leptonic acceleration mecha- 
nism, to distinguish it from a different type of acceleration, the hadronic one, to be 
discussed below, and which is believed by many as being the main mechanism for 
extragalactic ultra-high energy cosmic rays. Notably, AGNs are also extragalactic 
objects and hence hadronic acceleration mechanism may be relevant for the pro- 
duction of very high energy Gamma rays, as alternative scenarios to (or co-existing 
with ) SSC-lcptonic acceleration mechanism described above. 
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Fig. 11. A hadronic acceleration model for the production of very high energy gamma rays from 
cxtragalactic sources. High-energy Protons which have been accelerated in AGNs or other sources 
in the presence of magnetic cosmic fields, can interact with the protons of the interstellar medium 
to give rise to pions (neutral ir° and charged n^). Photons from neutral pion decays(7r° — y 27) 
could be detected together with those coming from charged-pion-to-muon conversion processes 
— y fjb — y e , which eventually yield electrons or positrons that scatter a la Inverse Compton 
with low-energy CMB photons resulting in detectable (on Earth or satellite experiments) high- 
energy photons (picture taken from H. Tajima SLAC-DOE (USA) Programme review talk (June 
7, 2006), (www-conf.slac.stanford.edu/programreview/2006/Talks) ). 



Hadronic Acceleration: A prominent way of producing high energy Gamma Rays 
of extragalactic origin, is the scattering of very high energy protons (produced in 
the jet of the AGNs by means of gravitational energy conversion mentioned above 
(c.f. fig-9)) off protons in the interstellar medium. Such collisions result in pion 
production, of which neutral pions decay (ir° — > 27) and give rise to very high 
energy photons that are detected directly. The charged pions on the other hand are 
converted to muons, whose decays produces electrons or positrons (% —y (1 —y 
e ); the scattering of the latter with low-energy photons of the cosmic microwave 
background radiation (CMB) then results, through (inverse) Compton scattering, 
in photons being detected on Earth or satellites (c.f. fig. 11). 
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Fig. 12. Comparing the spectra of very high energy gamma rays in models of leptonic and 
hadronic acceleration of cosmic particles. In the figure, E denotes the photon energy (in some 
generic units) and N the observed cosmic photon number, while the suffix "e" ("p") denotes quan- 
tities pertaining to electrons (protons) and B denotes the magnetic field. From the differences in 
the shape and position of these spectra one can get information on the kind of acceleration that 
takes place (picture taken from talk of M. de Naurois, Workshop of HSSHEP, April 2008, Olympia 
(Greece), (http : //www. inp . demokritos .gr /confer ences/HEP2008-01ympia/)). 



The energy spectra of hadronic acceleration models are different in shape and 
location of their peaks on the energy axis from those of leptonic acceleration, as can 
be seen in fig. 12, where various photon spectra in the Universe are superimposed 
for comparison. It is worthy of mentioning that in order to interpret the current high 
energy gamma ray data using leptonic inverse Compton measurements, one has to 
assume that the galactic magnetic fields are of low intensity. From measurements 
of high energy Gamma ray spectra from AGNs or other extragalactic sources, such 
as Gamma Ray Bursts (GRB), we can then soon get sufficient information into the 
precise way of acceleration of cosmic particles. The current experimental knowledge 
on cosmic high-energy gamma ray spectra can be summarised in fig. 13, from which 
it is clear that we need more measurements in the lower-energy part of the spectrum 
before conclusions can be drawn on the kind of cosmic acceleration taking place at 
various celestial sources. 
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Fig. 13. Recent data on Gamma ray spectra. The symbols in the axes have the same mean- 
ing as in fig. 12. To be sure on the kind of cosmic acceleration taking place (i.e. leptonic 
or hadronic) one needs data at the lower energy part of the spectrum. Such data can be 
provided, for instance, by the FERMI (GLAST) satellite, which was recently launched (pic- 
ture taken from talk of M. de Naurois, Workshop of HSSHEP, April 2008, Olympia (Greece), 
(http : //www. inp . demokritos.gr/conferences/HEP2008-01ympia/)). 

4.2. The MAGIC delays: Adding to the uncertainties on the VHE 
Gamma-Ray production mechanisms 

The observed time delays of order 4 ± 1 minutes between the most energetic (TeV) 
Gamma rays from AGN Mkn 501 (c./. fig. 3), observed by the MAGIC Telescope 3 ' 12 , 
lead to further uncertainties in the production mechanism of such photons. 

As discussed in 3 , the conventional model of SSC used to explain the origin of 
VHE Gamma Rays as being due to an electronic uniform acceleration in the AGN 
jet region, which finds a good application in other AGNs, such as Crab Nebula, 
fails to account for the observed time delay by MAGIC. The use of the acceleration 
parameters in the Crab Nebula AGN leads to only millisecond delays of the more 
energetic photons if applied to the Mkn 501 case. 

This prompted speculations that the conventional SSC mechanisms involving 
uniform acceleration of the relativistic blob of particles in the jet of the AGN (c.f. 
fig. 9) need to be modified for the Mkn 501 case. In fact, several propositions along 
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this line have been made so far: 

• (i) Particles inside the emission region moving with constant Doppler factor 
need some time to be accelerated to energies that enable them to produce 
7 rays with specific high energies in the TeV region 3 . 

• (ii) The 7-ray emission has been captured in the initial phase of the accel- 
eration of the relativistic material blob in the jet of the AGN (c.f. fig. 9), 
which at any point in time radiates up to highest 7-ray energies possible 82 . 

• (iii) A version of the SSC scenario (termed one-zone SSC model), which 
invokes a brief episode of increased particle injection at low energies 83 . 
Subsequently, the particles are accelerated to high energies, which thus 
accounts for the observed delays, but they also emit synchrotron and SSC, 
thereby loosing energy. As in scenario (ii) above, also according to this 
model the MAGIC observations have caught the relativistic electrons in 
the jet of the AGN at their acceleration phase. 

To the above I would also like to add the possibilities that some hadronic mech- 
anisms might also be in operation here, which seems not to have been discussed by 
the community so far. It therefore becomes clear from the above brief discussion 
that the situation concerning the delayed production of VHE Gamma Rays from the 
AGN Mkn 501 is far from being resolved by means of conventional (astro)physics 
at the source. 

This brings us to the main topic of this article, which is a possible (albeit specu- 
lative at this stage) link of the MAGIC observation with more fundamental physics 
associated with the very structure of space time on which the propagation of the 
VHE Gamma Rays takes place. 



4.3. Quantum- Gravity Space-Time Foam and the MAGIC delays: 
wild speculation or realistic possibility? 

In 12 it has been observed that, as a result of the ability of the experiment to mea- 
sure individual (within the accuracy of the observations of course) photons from 
Mkn 501, of various energies, it should be possible to reconstruct the peak of the 
flare of July 9th 2005 using dispersion relations of these individual photons during 
their journey from emission till observation. In fact, we went one step further and 
assumed sub-luminal modified dispersion relations of the type expected 33 > 41 > 13 > 15 
to be encountered in a model of quantum-gravity (QG) induced space-time foam 1 
coming from string theory 18 . As we have already discussed in the previous section, 
the sub-luminality of the QG-induced refractive index in such models is guaran- 
teed by the very nature of string theory, which respects the cornerstone of special 
relativity that the speed of light in vacuo is the maximal material velocity. In this 
respect the space-time foam in such theories leads to the absence of birefringence, 
in other words the refractive index is the same for both photon polarizations. This 
is an important feature, which allows the MAGIC results to be compatible with 
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other stringent limits of Lorentz invariance from other astrophysical sources, as we 
shall explain below. 

Although, as we have already argued, the observed time delays within the string 
model pertain to stringy uncertainties 18 , and as such may be disentangled from 
modified dispersion relations due to a breakdown of the local effective lagrangian 
formalism, nevertheless it is instructive to follow a historical path, and discuss first 
constraints on the photon dispersion relations coming from best fits to the observed 
delays and then proceed to the realistic string theory case. This will help the reader 
understand the failure of the former approach and the advantage of the string model 
over many approaches to QG based on local effective field theories. 

In 12 we examined two cases of QG-induced modified dispersion for photons, 
stemming from (16) upon assuming: 

• Case I: Photon Refractive index suppressed Linearly by the QG energy 
scale Mqqi, i.e. only the coefficient c\ > in the series of Eq. (16) is non 
zero, its positivity being required by the sub-luminal nature assumed for 
the propagation, as ensured by the string theory underlying model 33 > 15 > 18 . 

• Case II: Photon Refractive index suppressed Quadratically by the QG en- 
ergy scale Mqq2, *-e. only the coefficient c-i > in the series of Eq. (16) 
is non zero, its positivity again being linked to the sub-luminal nature as- 
sumed for the propagation. String theory models with this kind of quadratic 
suppression also exist in the modern approach to string theory, including 
representation of our world as a brane (domain wall hyperplane) 84 , but 
will not be discussed here. 

The method of reconstructing the peak ("most active part") of the flare by 
implementing modified dispersion relations for individual photons is based on the 
following well known fact of classical electrodynamics 85 : a pulse of electromagnetic 
radiation propagating through a linearly-dispersive medium, as postulated above, 
becomes diluted so that its power (the energy per unit time) decreases. The appli- 
cability of classical electrodynamics for estimating the low-energy behavior induced 
by space-time foam and the corresponding pulse-broadening effect have been dis- 
cussed in ref. 15 , where we refer the interested reader for further details and explicit 
examples. 

The dilution effects for the linear or quadratic cases may easily be obtained as 
described in 85 by applying the dispersion laws 

Lo(k) = fc[l - fc/(2M QG i) , or u(k) = k[l - k 2 /(3M^ G2 )} , (82) 

where to denotes the frequency of the photon, with wave vector k. Any transforma- 
tion of a signal to reproduce the undispersed signal tends to recover the original 
power of the pulse. If the parameter Mqgii , n = 1, 2 is chosen correctly, the power 
of the recovered pulse is maximized. 

This observation has been implemented in the analysis of 12 by appropriately 
choosing (using statistical-analysis techniques) a time interval (ti; t2) containing the 
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most active part of the flare. For the record, we mention that this procedure has 
been applied in 3 to 1000 Monte Carlo (MC) data samples generated by applying 
to the measured photon energies the (energy-dependent) Gaussian measurement 
errors. The results of the reconstruction of the peak ("most active part") of the 
flare using the linear Case I are demonstrated in fig. 14 for completeness. In a 
similar way one gets bounds on the quadratic Case II. These results have also been 
confirmed using different statistical analysis techniques, independent of the ECF. 
The interested reader is referred to 12 for details of the analysis. 
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Fig. 14. The Left figure shows the Energy Cost Functional (ECF) from one realization of the 
MAGIC measurements with photon energies smeared by Monte Carlo, for the case of a vacuum 
refractive index that is linear in the photon energy (Case I). The ECF exhibits a clear maximum, 
whose position may be estimated by fitting it with a Gaussian profile in the peak vicinity. The 
Right figure shows the results of such fits to the ECFs, specifically the t; distribution for the 1000 
energy-smeared realizations of the July 9 flare, where t; is defined through the representation 
of the QG scale as Mqqi = 1.445 X 10 16 s/t;. From this distribution we derive the value t; = 
(0.030 ± 0.012) s/GeV, leading to a lower limit M QG1 > 0.21 X 10 18 GcV at the 95% C.L. 



Taking into account the uncertainties in the production mechanism at the source, 
which could also contribute as we have discussed before to the observed delays, we 
can only place lower bounds on the quantum gravity energy scale from such an 
analysis. In fact for the linear and quadratic cases we obtain the lower bounds (c.f. 
fig 14): 

Mqqi > 0.21 x 10 18 GeV at 95% Confidence Level (C.L.), 

M QG2 > 0.26 x 10 11 GeV at 95% C.L.. (83) 

It is important to notice that, had the mechanism at the source been understood, 
and the emission of the different photons been more or less simultaneous, i.e. by at 
least two orders of magnitude smaller as compared to the observed delays, then the 
above lower bounds could be turned into a real measurement of the quantum gravity 
scale. It is surprising, therefore, that in such a situation Case I can reproduce the 
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observed delays, provided Mqqi is of order of the so-called reduced Planck mass, 
which is an energy scale characterizing conventional string theory models. 

In the analysis of 12 it was possible to exclude the possibility that the observed 
time delay may be due to a conventional QED plasma refraction effect induced as 
photons propagate through the source. From the discussion in sub-section 2.1, in 
particular eq. (11), it becomes clear that if the delay would be due to plasma effects 
at the source region, then this would induce 

At = D(a 2 T 2 /6k 2 ) In 2 (kT/m 2 e ) , 

where a is the fine-structure constant, k is the photon momentum, T is the plasma 
temperature, m e is the mass of electron, D is the size of the plasma, and we use 
natural units: c,h~ 1. Plausible numbers such as T <~ 1CP 2 MeV and D <~ 10 9 km 
(for a review see 86 ) yield a negligible effect for k <~ 1 TeV, which is of order 
of the photon momenta relevant to the MAGIC experiment. Exclusion of other 
source effects, such as time evolution in the mean emitted photon energy, might be 
possible with the observation of more flares, e.g., of different AGNs at varying 
redshifts. Observations of a single flare cannot distinguish the quantum-gravity 
scenarios considered here from modified synchrotron-self-Compton mechanisms. 

4.4. H.E.S.S. and FERMI Observations and Quantum- Gravity 
Scale Bounds 

However, the above-described pioneering study demonstrates clearly the potential 
scientific value of an analysis of multiple flares from different sources. The most 
promising candidate for applying the analyses proposed here is the flare from 
the Active Galaxy PKS 2155-304 detected recently by H.E.S.S. (tfigh Energy 
Stereoscopic System) Collaboration 9 , another experiment involving arrays of 
Cherenkov Telescopes in Namibia (Africa). This galaxy lies further than Mk501 
at redshift z = 0.116 and there is a much higher statistics of photons at energy 
ranges of a few TeV. 

In fact H.E.S.S. collaboration published their measurements 10 on the arrival 
time of photons from PKS 2155-304. However, unlike the MAGIC observations 
from Mk501 Galaxy, there was no time lag found between higher- and lower-energy 
photons in this case. These results can thus place only bounds on the quantum 
gravity scale, if space-time foam is assumed to affect the photon propagation. The 
bounds are similar to the MAGIC case (83) above. 

The H.E.S.S. result can mean several things, and certainly points towards differ- 
ent source mechanisms for the acceleration of photons between the two galaxies Mkn 
501 and PKS 2155-304. However, this second measurement by H.E.S.S. cannot still 
rule out the possibility that quantum gravity plays a role in photon propagation. 
For instance, one cannot exclude the (admittedly remote) possibility that, due to a 
still unknown source effect, the high-energy photons in the PKS 2155-304 Galaxy 
are emitted first, in contrast to the Mkn 501 case, in such a way that a sub-luminal 
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vacuum refractive index quantum-gravity effect, of a strength appropriate to pro- 
duce the delays observed by the MAGIC experiment, "slowed these photons down" 
as compared to the lower-energy ones, so that there are no observable delays in the 
arrival times between the higher and lower energy photons in the H.E.S.S. experi- 
ment. Of course, one cannot exclude the possibility that the conditions of this set 
of measurements, for some reason, prohibited the detection of an observable time 
lag between high and low energy photons. 



GRB 09092B Fermi 09/2009 
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Fig. 15. Comparison of data on delays At in the the arrival times of energetic gamma rays from 
various astrophysical sources with models in which the velocity of light is reduced by an amount 
linear in the photon energy 8 . The graph plots on a logarithmic scale the guantity At/E and a 
function of the red-shift, K(z), which is essentially the distance of the source from the observation 
point. The data include two AGNs, Mkn 501 3,12 and PKS 2155-304 10 , and three GRBs observed 
by the Fermi satellite, 090510 6 , 09092B 7 and 080916c 5 . 

Although the available sample of AGNs is still not large enough for a robust 
analysis on bounds of Quantum-Gravity medium effects 13 > 15 , nevertheless, one can 
at least check for consistency between the available MAGIC and HESS results, 
and gauge the magnitude of possible intrinsic fluctuations in the AGN time-lags. 
Comparing the time-lag measured by MAGIC for Mkn 501 at redshift z = 0.034: 
At/E y = 0.030 ±0.012 s/GeV, with that measured for PKS 2155-304 at z = 0.116: 
At/E 7 — 0.030 ± 0.027 s/GeV, we see that they are compatible with a common, 
energy-dependent intrinsic time-lag at the source. On the other hand, they are also 
compatible with a universal redshift- and energy-dependent propagation effect: 

At/R r = (0.43 ± 0.19) xtf(*)B/GeV, K{z) = f ^% z)3 , (84) 
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assuming an expanding Universe within the framework of the standard 
Cosmological-constant-Cold-Dark-Matter (ACDM) model. The best fit 8 of the 
MAGIC and HESS data based on (84) leads to the following result for the 
Quantum Gravity scale, assuming that it is the dominant cause of the delay: 
M QG1 = (0.98±g;|5) x 10 18 GeV. The situation is depicted in fig. 15. 

With measurements from only a few available flares from AGN it is, therefore, 
not possible to disentangle with any certainty source from propagation effects. For 
this we need statistically significant populations of available data 15 . Unfortunately 
the occurrence of fast flares in AGNs is currently unpredictable, and since no corre- 
lation has yet been established with observations in other energy bands that could 
be used as a trigger signal, only serendipitous detections are currently possible. It 
seems unlikely that the relatively rare and unpredictable sharp energetic flares pro- 
duced only occasionally by AGNs, which have a relatively restricted redshift range 
and hence a small lever arm, will soon be able to provide the desired discrimination. 

On the other hand, Gamma Ray Bursts (GRBs) are observed at a relatively 
high rate, about one a day, and generally have considerably larger redshifts. The 
advent of the FERMI (nee GLAST) Telescope with its large acceptance offers the 
possibility of achieving the required sensitivity. Indeed, the FERMI Collaboration 
has already made a report of GeV-range 7 rays from the GRB 080916c 5 . In this 
GRB, there is a 4.5-second time-lag between the onsets of high- (> 100 McV) 
and low-energy (< 100 KeV) emissions. Moreover, the highest-energy photon GRB 
080916c measured by the FERMI 7-ray telescope had an energy E = 13.2 till GeV, 
and was detected At = 16.5 s after the start of the burst. Spectroscopic information 
has been used by the GROND Collaboration 87 to estimate the redshift of GRB 
080916c as z = 4.2 ± 0.3 5 . Assuming this value of the redshift, the best fit (84) 
would correspond to a time-lag 

At = 25 ills (85) 

for a 13 GeV photon from GRB 080916c. As discussed in 8,5 , such time de- 
lays can fit excellently within the above-mentioned QG scenario for a sublumi- 
nal refractive index for photons, with the following lower bound for the QG scale 
Mqgi > 1-50 ± 0.20 x 10 18 GeV, where the inequality is due to the ignorance 
of the source mechanism (c.f. fig. 15). This bound is consistent with the MAGIC 
and HESS results stated previously. The reader should also bear in mind that the 
4.5-second time-lag observed for ~ 100 MeV photons could not be explained by a 
propagation effect that depends linearly on the energy 5 . It is clear therefore that 
the FERMI Telescope has already demonstrated the sensitivity to probe a possible 
linearly energy- dependent propagation effect at the level reached by the available 
AGN data, and it is appropriate and possibly helpful to consider how such an effect 
could be probed in the future. 

The analysis of 8 has also demonstrated that although the three above-mentioned 
sets of data, from MAGIC, FERMI and HESS Collaborations, can be explained 
simultaneously by a linear in energy vacuum refractive index, suppressed by a single 
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power of the quantum gravity scale, this is not the case for a refractive index scaling 
quadratically with the photon energy, for instance of the type encountered in some 
brane models with asymmetric warp factors, as in ref. 84 . Indeed, on assuming that 
the quadratic refractive index is the sole cause of the observed delays in the arrival 
of high- vs. low- energy photons in both the MAGIC and FERMI cases, this would 
imply a time delay of order 0.24± 0.16 s for the most energetic photon (13.22 GeV) 
of the GRB 080916c, for a quantum gravity scale that saturates Mqg2 in (83). 
This is two orders of magnitude smaller than the measured time-lag (16.5 s) by the 
FERMI Collaboration 5 . 

4.5. Trouble with GRB 090510 - non-uniform D-particle foam? 

In figure 15 we also depict the data from the short burst GRB 090510 6 , observed 
by the FERMI Telescope, which so far we did not include in our linear in energy 
fit. It is clear that a fit with a linearly modified dispersion relation does not work 
in this case for the values of the Quantum Gravity scale Mqq f lc ~ 10 18 GeV that 
fit the MAGIC, H.E.S.S. and the other FERMI data. Indeed, the observed short 
delays of this burst can be explained on the basis of linearly modified dispersion 
relations only if the quantum gravity scale is larger than Mqgi — 1.2Mp where 
Mp = 1.22 x 10 19 GeV is the Planck mass, which thus is incompatible with the value 
of the quantum gravity scale for the MAGIC data fit. This has been used by the 
Fermi collaboration as a strong indication that Lorentz violating models entailing 
modified dispersion relations with linear suppression are ruled out on naturalness 
grounds, since the so-obtained quantum gravity scale for the GRB 090510 is about 
1.2 times larger than the four-dimensional Planck mass Mp. 

We would not agree with this statement. Leaving aside the fact that from a single 
measurement, with uncertainties on the precursor of the GRB, one cannot draw safe 
conclusions, we mention that what one calls a natural scale of QG is highly model 
dependent. For instance, as we have seen in the previous subsection, in the string 
foam model, the relevant "QG scale" , that dictates the order of the foam-induced 
time delays of photons, is a complex function of the model's parameters and is not 
simply given by the string scale M s 11 . Indeed, as becomes clear from (79) (or (81)) 
the relevant scale is not simply the D-particle mass, M s /g s , but a combination 

M Q G-D-foam ~ 9e n'\z) ( 86 ) 

involving the linear number density n*(z) of the foam defect, encountered by the 
photon during its propagation from the source till observation. This quantity de- 
pends on the bulk density of the D-particles, which in the model of 17 is a free 
parameter. 

Inhomogeneous bulks are perfectly consistent background configurations for our 
brane world scenario 11 . In order to match the photon delay data of Fig. 15 with 
the D-foam model, we need a reduction of the linear density of defects encountered 
by the photon by about two orders of magnitude in the region 0.2 < z < 1, whereas 



String Quantum Gravity, LIV and Gamma-Ray Astronomy 63 

for z < 0.2 there must be, on average, one D-particle defect per unit string length 
l s . We therefore assume that our D-brane encountered a D-void when 0.2 < z < 1, 
in which there was a significant reduction in the bulk nine-dimensional density of 
defects. Such assumptions can be consistent with cosmological considerations on the 
dark sector of the model, as explained in 11 . In particular, an important parameter 
for the cosmology of the model is the propagation velocity v of the brane world 
in the bulk (c.f. fig. 6). As we have discussed previously (48) this is related 11 to 
the density of the defects near the brane world, n short (z), and hence to the linear 
density of defects encountered by the photon on the brane world. 

To have a reduced density by two order of magnitude at redshifts z = 0.9, 
while having a density of defects of 0(1) per string length at redshifts z < 0.1, 
one must consider the magnitude of v as well as the string scale t s . For instance, 
it follows from (48) that for string energy scales of the order of TeV, i.e. string 
time scales l s /c = 10 -27 s, one must consider a brane velocity v < VTO x 10~ n c, 
which is not implausible for a slowly moving D-brane at a late era of the Universe 
This is compatible with the constraint on v obtained from inflation in 57 , namely 
v 2 < 1.48 x lO -5 .^ 1 , where g s < 1 for the weak string coupling we assume here. In 
our model, due to the friction induced on the D-brane by the bulk D-particles, one 
would expect that the late-epoch brane velocity should be much smaller than that 
during the inflationary era immediately following a D-brane collision 57 . 

The above considerations can be extended 11 to type-IIB constructions of D- 
foam 19 . If there was a depletion of D-particle defects in a certain range of redshifts 
in the past, e.g., when z <~ 0.9, the volume Va3 in (66) would have been larger than 
the corresponding volume at z <C 1, and the corresponding coupling reduced, and 
hence also the corresponding cross section a describing the probability of interaction 
of a photon with the D-particle in the foam. 

The effective linear density of defects n*{z) appearing in (79) can then be related 
to the number density per unit three- volume on the compactificd D7-brane via 
n* (z) ~ any'. The quantity is, in turn, directly related to the bulk density 
of D-particles in this model, where there is no capture of defects by the moving 
D-brane, as a result of the repulsive forces between them. The cross section a is 
proportional to the square of the string amplitude (70) describing the capture of an 
open photon string state by the D-particle in the model, and hence is proportional to 
g| 7 , where g^7 is given in (66). If recoil of the D-particle is included, the amplitude is 
proportional to the effective string coupling: g^ s oc gf 7 , where g eS = g s (l — \u\ 2 y^ 2 , 
with u <C 1 the recoil velocity of the D-particle 18 > 19 . One may normalize the string 
couplings in such a way that n*(z = 0) = 0(1), in which case n*(z ~ 0.9) <C 1 
by about two orders of magnitude. However, in general the induced gauge string 
couplings depend crucially on the details of the compactification as well as the 
Standard Model phenomenology on the D-brane world. Hence the precise magnitude 



'Much smaller velocities are required for small string scales that are comparable to the four- 
dimensional Planck length. 
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of the time delays is model-dependent, and can differ between models. 

The above discussion hopefully demonstrates that the question: "what is a nat- 
ural scale for QG?" is highly dependent on the details of the microscopic model 
and does not admit a straightforward answer. Many more measurements, at various 
cosmological distances (red shifts) are needed before we falsify a model, such as the 
string foam model, for instance, discussed here. Nevertheless, with the progress in 
extraterrestrial instrumentation, and with the launch of facilities like FERMI or 
future Cherenkov Telescope Arrays (CTA) 88 , this will hopefully become possible 
in the foreseeable future. 

4.6. Other (astrophysical) constraints on quantum- gravity foam 

The sensitivity of the MAGIC (and FERMI) observations to Planck scale physics 
(83), at least for linearly suppressed modified dispersion relations, calls for an im- 
mediate comparison with other sensitive probes of non-trivial optical properties of 
QG medium. 

Indeed, from the analysis of 12 , there was no microscopic model dependence 
of the induced modifications of the photon dispersion relations, other than the 
sub-luminal nature of the induced refractive index and the associated absence of 
birefringence, that is the independence of the refractive index on the photon po- 
larization. The latter feature avoids the otherwise very stringent constraints on the 
photon dispersion relation imposed by astrophysical observations, as we now come 
to discuss. 

We shall be very brief in our description of the complementary astrophysical 
tests on Lorentz invariance and quantum-gravity modified dispersion relations, to 
avoid large diversion from our main point of this review article which is string 
theory. 

There are three major classes of complementary astrophysical constraints, to 
be considered in any attempt to interpret the MAGIC, FERMI or more general 7- 
ray Astrophysics results in terms of quantum-gravity induced anomalies in photon 
dispersion. 

• Birefringence and strong constraints on QG-induced photon dispersion: 
In certain models of quantum gravity, with modified dispersion relations, 
for instance the so-called loop-quantum gravity 36 , the ground state breaks 
reflexion symmetry (parity) and this is one of the pre-requisites for a de- 
pendence of the induced refractive index on the photon polarization, i.e. 
birefringence. We remind the reader that in birefringent materials this is 
caused precisely by the existence of some kind of anisotropics in the mate- 
rial. The velocities of the two photon polarizations (denoted by ±) in such 
QG models may be parametrised by: 

v± =c(l±£(?vuj/Mp) n ) (87) 

where Mp = 1.22 x 10 19 GeV is the Planck energy scale, and £ is a parameter 
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following from the underlying theoretical model, which is related with the 
modifications of the pertinent dispersion relations for photons. The order 
of suppression of these effects is described by n which in the models of 
36 assumed the value n = 1, but in general one could have higher order 
suppression, as we have discussed in (16). 

Vacuum QG birefringence should have showed up in optical measurements 
from remote astrophysical sources, in particular GRBs. Ultraviolet (UV) ra- 
diation measurements from distant galaxies 71 and UV/optical polarization 
measurements of light from Gamma Ray Bursters 72 rule out birefringence 
unless it is induced at a scale (way) beyond the Planck mass (for linear mod- 
els, the lower bound on the QG scale in such models can exceed the Planck 
scale (~ 10 19 GeV) by as much as seven orders of magnitude). Indeed, in 
terms of the parameter n introduced above (c.f. (87)), for the case n = 1 
of 36 , one finds from optical polarization observations that the absence of 
detectable birefringence effects imply the upper bound |£| < 2 x 10~ 7 , which 
is incompatible with the MAGIC observed delays, saturating from below 
the bounds (83). 

At this point, we wish to mention that, using recent polarimetric observa- 
tions of the Crab Nebula in the hard X-ray band by INTEGRAL 89 , the 
authors of 73 have demonstrated that the absence of vacuum birefringence 
effects constrains linearly suppressed Lorentz violation in quantum electro- 
dynamics to the level |£| < 6 x 1CP 10 at 95% C.L., thereby tightening by 
about three orders of magnitude the above-mentioned constraint. 
• Synchrotron Radiation and further stringent constraints for electronic QG- 
induced anomalous dispersion in vacuo. 





Fig. 16. The Crab Nebula (left image) is a supernova remnant, with a rotating neu- 
tron star (the Crab Pulsar) at its centre (middle image). Observations of synchrotron ra- 
diation (right image, radiation spectrum (arbitrary units)) from such celestial objects place 
very stringent constraints on quantum gravity models with anomalous dispersion relations for 
electrons (images taken from NASA, http://www.nasa.gov/ (Crab Nebula, left), wikipedia, 
http : //en . wikipedia . org/wiki/CrabNebula (Crab Pulsar, middle) and the first of Ref. 45 (right)) . 



Another important experimental constraint on models with QG-induced 
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anomalous dispersion relations comes from observations of synchrotron ra- 
diation from distant galaxies 74 ' 45 > 75 7 such as Crab Nebula (c.f. fig. 16). As 
mentioned previously, the magnetic fields at the core regions of galaxies 
curve the paths of (and thus accelerate) charged particles, in particularly 
electrons (which are stable and therefore appropriate for astrophysical ob- 
servations), and thus, on account of energy conservation this results in 
synchrotron radiation. 

In standard electrodynamics 85 , electrons in an external magnetic field H, 
follow helical orbits transverse to the direction of H. The so- accelerated 
electrons in a magnetic field emit synchrotron radiation with a spectrum 
that cuts off sharply at a frequency uj c (c.f. fig. 16, right panel): 

, ,LI _ 3 eH 1 (oo\ 

where e is the electron charge, mo its mass, and f3± = v± is the component 
of the velocity of the electron perpendicular to the direction of the magnetic 
field. The superfix LI in (88) stresses that this formula is based on a LI 
approach, in which one calculates the electron trajectory in a given magnetic 
field H and the radiation produced by a given current, using the relativistic 
relation between energy and velocity. 

All these assumptions are affected by violations of Lorentz symmetry, such 
as those encountered in quantum-gravity space-time foam models, leading 
to modified dispersion relations of the form: 

W 2 (fc) = fc 2 +£ 7 ^, (89) 

E 2 (p)=ml+p 2 +^^-, (90) 

for photons (89) and electrons (90), where uj and k are the photon frequency 
and wave number, and E and p are the electron energy and momentum, 
with mo the electron (rest) mass. In the spirit of the MAGIC observation 
analysis above, we assume here linear (a = 1) or quadratic QG (a = 2) 
effects, characterized by parameters £ 7 and £ e , extracting the Planck mass 
scale Mp = 1.22 x 10 19 GeV. In fact one can do the analysis 45 for a 
general a (single power) and attempt to extract limits on this parameter 
by matching with observations. 

A detailed analysis 74 > 45 J including the modifications in the electron's tra- 
jectories due to space-time foam 45 , yields: 

d + Vi^7v) 1 / 2 ^ + ( a + i)(^) c 

M = Mp/\^ e \, 77 = 1 - [E/MY , (91) 

where co^ 1 is given in (88) and the superscript "QG" indicates that the 
QG-modified dispersion relations (90) are used. This function is plotted 
schematically (for a = 1) in fig. 16 (right panel). 
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In 74 > 45 7 the above QG- modified dispersion relations have been tested using 
observations from Crab Nebula. It should be emphasized that the estimate 
of the end-point energy of the Crab synchrotron spectrum and of the mag- 
netic field used above are indirect values based on the predictions of the 
Synchrotron Self-Compton (SSC) model of very-high-energy emission from 
Crab Nebula 90 . In 45 the choice of parameters used was the one that gives 
good agreement between the experimental data on high-energy emission 
and the predictions of the SSC model 90 ' 86 . Estimating the magnetic field 
of Crab Nebula in the region 160 x 1CT 6 Gauss < H < 260 x 10~ 6 Gauss, 
and requiring |£ e | < 1 (which thus sets the quantum gravity scale as at least 
Mp) one obtains the following bounds on the exponent a of the dispersion 
relations (90) 45 : 

a > a c : 1.72 < a c < 1.74 (92) 

These results imply already a sensitivity to quadratic QG corrections with 
Planck mass suppression Mp. 

However, for photons there are no strong constraints on £ 7 coming from 
synchrotron radiation studies, unless in cases where QG models entail bire- 
fringence 75 , where, as we discussed above, strong constraints on the photon 
dispersion are expected at any rate from optical measurements on GRBs. 
In this sense, the result (92) excludes the possibility that the MAGIC ob- 
servations leading to a four-minute delay of the most energetic photons are 
due to a quantum foam that acts universally among photons and electrons. 
However, the synchrotron radiation measurements cannot exclude anoma- 
lous photon dispersion with linear Planck-mass suppression, leading to a 
saturation of the lower bound (83), in models where the foam is transpar- 
ent to electrons, as in the string foam case. 
• Strong constraints from Ultra-high- energy Cosmic photon annihilation 
Further strong constraints on generic modified dispersion relations for pho- 
tons, like the ones used in the aforementioned QG-interpretation of the 
MAGIC results 12 , comes from processes of scattering of ultra- high-energy 
photons, with energies above 10 19 eV, off very-low energy cosmic photons, 
such as the ones of the cosmic microwave background (CMB) radiation that 
populates the Universe today, as a remnant from the Big-Bang epoch. In 
20 it has been argued that the non-observation of such ultra-high energy 
(UHE) photons places very strong constraints on the parameters governing 
the modification of the photon dispersion relations, that are several orders 
of magnitude smaller than the values required to reproduce the MAGIC 
time delays, should the effect be attributed predominantly to photon prop- 
agation in a QG dispersive medium. 

The main argument relies on the fact that an ultra-high-energy photon 
would interact with a low-energy ("infrared") photon of the CMB back- 
ground (with energies in the eV range) to produce electron prositron pairs, 
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according to the reaction: 

7uhe + 7cmb => e + er . (93) 

The basic assumption in the analysis is the strict energy and momentum 
conservation in the above reaction, despite the modified dispersion rela- 
tions for the photons. Such an assumption stems from the validity of a 
local-effective-lagrangian description of QG foam effects on particles with 
energies much lower than the QG energy scale (assumed close to Planck 
scale Mp = 10 19 GeV). In this formalism, one can represent effectively the 
foam dispersive effects by higher-derivative local operators in a flat-space- 
time Lagrangian. The upshot of this is the modification of the pertinent 
equations of motion for the photon field (which in a Lorentz-invariant the- 
ory would be the ordinary Maxwell equations) by higher-derivative terms, 
suppressed by some power of the QG mass scale. 

One considers the modified dispersion relations (89), (90), which in the 
notation of 20 , taking explicit account of the various polarizations and he- 
licities, can be written as: 

El ± =pl + ml+rt ± pl(^) (94) 

with (u>, k) the four-momenta for photons, and (E,p e the corresponding 
four-momentum vectors for electrons; the suffix b indicates a low-energy 
CMB photon, whose dispersion relations are assumed approximately the 
normal ones, as any QG correction is negligible due to the low values of 
energy and momenta. The +(-) signs indicate left(right) polarizations (pho- 
tons) or hclicities (electrons). Positive (negative) £ indicate subluminal (su- 
perluminal) refractive indices. Upon the assumption of energy-momentum 
conservation in the process, one arrives at kinematic equations for the 
threshold of the reaction (93), that is the minimum energy of the high- 
energy photon required to produce the electron-positron pairs. 
For the linear- or quadratic- suppression case (Cases I and II in the MAGIC 
analysis above, for which n = 1, 2 respectively (94)), one finds that, for the 
relevant subluminal photon refractive indices corresponding to the satura- 
tion of the Mqoi lower bound in (83), the threshold for pair production dis- 
appears for ultra-high-energy photons, and hence such photons should have 
been observed. The non-observation of such photons implies constraints for 
the relevant parameters £, n which are stronger by several orders of magni- 
tude than the bounds (83) inferred from the MAGIC observations. 
From the analysis of 20 one concludes that in the case of linear Planck-mass 
suppression of the sub-luminal QG-induced modified dispersion relations 
for photons, of interest for the QG-foam interpretation of the MAGIC re- 
sults 12 , parameters with size £i > 10 14 are ruled out. This exceeds the 
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sensitivity of the MAGIC experiment (c.f. (83)) to such Lorentz-symmetry 
violating effects by fifteen orders of magnitude ! Similar strong constraints 
are also obtained from the non observations of photon decay (7 — > e + e~), 
a process which, if there are modified dispersion relations, is in general 
allowed 20 . 



4.7. The String Foam Models Evade the above Constraints 

From the above discussion it becomes clear that any model of refraction in space- 
time foam that exhibits effects at the level of the MAGIC sensitivity 12 (83) should 
be characterised by the following specific properties: 

• (i) photons are stable (i.e. do not decay) but should exhibit a modified 
subluminal refractive index with Lorentz-violating corrections that grow 
linearly with E/(Mqq 1 c 2 ), where Mqq 7 is close to the Planck scale, 

• (ii) the medium should not refract electrons, so as to avoid the synchrotron- 
radiation constraints 74 > 45 ; and 

• (iii) the coupling of the photons to the medium must be independent of 
photon polarization, so as not to have birefringence, thus avoiding the per- 
tinent stringent constraints 71 < 72 < 75 < 73 . 

• (iv) The formalism of local effective lagrangians should break down, in the 
sense that there are quantum fluctuations in the total energy in particle 
interactions, due to the presence of a quantum gravitational 'environment', 
such that stringent constraints, which otherwise would have been imposed 
from the non-observation of ultra-high energy photons (hu> > 10 19 eV), are 
evaded. 

The string- foam models (both type II A and type IIB) are characterized by all 
these properties 17 ' 45 > 18 ; and thus avoid the stringent constraints. The absence of 
birefringence and the transparency of foam to electrons or charged probes, for rea- 
sons of charge conservation, make the models surviving the stringent constraints 
from synchrotron radiation from distant Nebulae 74 ' 75 . 

Moreover, the analysis in 20 is based on exact energy momentum conservation 
in the process (93), stemming from the assumption of the local-effective lagrangian 
formalism for QG-foam. As we discussed in 91 > 22 ; and mentioned in section 3.7, 
however, such a formalism is not applicable in the case of the recoiling D-particle 
space-time foam model, where the fluctuations of space-time or other defects of 
gravitational nature paly the role of an external environment, resulting in energy 
fluctuations in the reaction (93). The presence of such fluctuations does affect the 
relevant energy-threshold equations, for the reaction to occur, which stem from 
kinematics, in such a way that the above stringent limits are no longer valid. In 
particular, we have seen that the pertinent anomalous dispersion terms are induced 
by the Finsler-like metric (58), which in turn arises by the distortion of space time 
due to the recoil of the D-particle space-time defect during its interaction with the 
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photon. This metric is quadratically suppressed by the string scale, and hence the so- 
induced modified dispersion relations (61) contain anomalous terms quadratically 
suppressed by the QG scale, a case for which the constraints coming from ultra 
high energy cosmic rays are weak. It should be stressed again that in the string 
model the modified dispersion relations are disentangled from the time delays (62) 
of the more energetic photons, which are linearly suppressed by the string scale. 
These delays are associated with the string uncertainty principle 49 and thus are 
not represented within the local effective field theory framework. We have also seen 
in section 4.5 that this disentanglement allows for a consistent interpretation of the 
MAGIC delays with those of the short GRB 090510 observed by FERMI within the 
string D-foam framework, provided of course the foam is inhomogeneous. 

In addition, as discussed in section 3.6, the form of the effective string coupling 
(59), and the induced metric (58), imply an upper bound in the photon momentum 
transfer (60), which characterizes the D-foam model as a result of the sub-luminal 
nature of the D-particle recoil velocity. These considerations imply that very high 
energy cosmic rays, when interacting with the D-foam, will have an extremely sup- 
pressed interaction rate (the string amplitudes are vanishing when the recoil velocity 
approaches the speed of light), thus providing additional reasons 22 for evading the 
strong constraints of 20 . There are also phenomenologically realistic brane models, 
with large extra dimensions 92 , for which the ratio M s /g s (i.e. in our context the 
D-particle mass) may be of the order of the conventional GZK cutoff 66 , 10 20 eV. In 
these models, in view of (60), there would be no photons with energies higher than 
this that could not be completely absorbed when interact with the D-particle foam, 
thereby providing an explanation for their absence, in accordance with observations. 

5. Conclusions and Outlook 

In this review we have discussed a stringy version of a (Lorentz-Invariance- Violating) 
space-time foam model as a candidate theory that provides an explanation for the 
delays of the more energetic photons from celestial sources, as observed by MAGIC 
and FERMI Telescopes, in agreement with all other current astrophysical tests of 
Lorentz Violation. This of course does not mean that there are no conventional 
astrophysics explanations for these delayed photon arrivals, but it demonstrates 
clearly that string theory (or better its modern version involving D-brane defects) 
is capable of explaining the observed photon delays, in agreement with all the other 
astrophysical data currently available. This is at least amusing, since it provides 
a framework for experimentally testing some models of string theory (entailing 
Lorentz Violation) at present or in the foreseeable future. 

The key point in the approach is the existence of space-time defects in the 
ground state of the model, whose topologically non-trivial interactions with the 
string states, via string-stretching during the capture process (c.f. fig. 7), are mainly 
responsible for the observed delays. The latter are found proportional to the inci- 
dent photon energy. The peculiarity of the D-particle foam in being transparent to 
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charged particles (as a result of electric charge conservation requirements), evades 
the stringent constraints on linear Planck scale suppression refractive indices that 
would be induced by electron synchrotron radiation studies from Crab Nebula 74 > 75 . 
Moreover, the absence of birefringence avoids the similarly stringent constraints on 
such models that otherwise would be imposed by galactic measurements 75 . Finally, 
it worths mentioning that the string-stretched linear in energy time delays (49), (81), 
when applied to neutrinos, can be flavour (i.e. neutrino-species) independent, thus 
avoiding 93 the stringent constraints that would be obtained from models of quan- 
tum gravity with flavour-dependent modifications of neutrino propagation and thus 
modifications in their oscillations 94 . 

As already mentioned, the MAGIC result needs confirmation by other experi- 
ments like H.E.S.S. 9 ' 10 , or other photon dedicated experiments, like FERMI (for- 
merly GLAST) 4 , where photons from Gamma Ray bursts will be observed. FERMI 
observed an extremely intense and short burst, GRB 090510, which imposes a strin- 
gent constraint on the allowed effective quantum gravity scale for linear anomalous 
refraction effects, Mqq > l.2Mp. Assuming that the measurement is correct, we 
could incorporate such a stringent constraint into our D-foam model, in a way 
consistent with the MAGIC and the other FERMI measurements on delayed pho- 
ton arrivals, by invoking inhomogeneous D-particle foam models, exhibiting voids 
around the red-shift region of the GRB 090510. Such models can provide consistent 
cosmologies within our framework. 

We must stress once more, however, that despite the apparent "success" of the 
D-foam model in dealing with stringent astrophysical constraints so far, one cannot 
as yet draw any safe conclusions on the validity of the model in Nature. Many more 
high-energy astrophysical photon measurements are needed in order to disentangle 
source from possible propagation effects due to fundamental physics. If a statis- 
tically significant population of data on photons from cosmic sources is collected, 
exhibiting refractive indices varying linearly with the distance of the source 15 , as 
well as the photon energy, then this would be a very strong confirmation of the 
D-particle foam model, for reasons explained above. However, it must be noted 
that GRB's, which are expected to lead to statistically significant data in the next 
few years, will produce photons much lower in energies than the flares observed in 
AGN, and this could be a drawback. At any rate there are attempts to claim that 
observations from FERMI will have sensitivity close to the Planck scale 95 for such 
linear-suppression models soon. The case of GRB090510 is a perfect example of how 
a single measurement of a short, intense high energy burst, can place stringent limits 
on Lorentz Violation at the Planck scale. However, caution should be exercised here 
as to what one means by "sensitivity at the Planck scale" . As demonstrated above 
(c.f. (79)), the effective quantum gravity scale is actually a complicated function of 
many microscopic parameters in the model. Nevertheless, if many observations at 
various rcdshift regimes on delayed arrivals of cosmic photons become available in 
the future, then we shall be able to make some definite conclusions regarding the 
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order of magnitude of possible quantum gravity effects and thus falsify models, such 
as the D-particle space-time foam. 

From the above discussion it becomes hopefully clear to the reader that experi- 
mental searches for quantum gravity, if the latter is viewed as a medium, are highly 
model dependent as far as the sensitivity to experimental falsification of the predic- 
tions of the underlying model is concerned. However, we are entering an era where 
low-energy (compared to Planck scale) physics experiments can already provide 
valuable information on the structure of space-time at the scales where Quantum 
Gravity is expected to set in. Very- and Ultra- high energy Astrophysics is at the 
forefront of such fundamental research. We therefore expect that, for the years to 
come, this branch of physics will proceed in parallel with terrestrial high energy ex- 
periments, such as the Large Hadron Collider launched at CERN recently, and be 
able to provide us soon with complementary important information on the under- 
lying fundamental structure of our Cosmos. Time will then show whether quantum 
gravity can be finally put to experimental confirmation. 
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